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Abstract 

We investigate 'non-local' schemes for baryogenesis at a first order electroweak phase 
transition, in which the effects of a CP violating condensate on the bubble wall propagate 
into the unbroken phase where the sphaleron rate is unsupressed. Such a condensate 
exists in multi-Higgs extensions of the standard model, and may exist due to an instability 
in the minimal standard model. In this paper we first discuss the general problem of 
determining the CP violating perturbations, distinguishing two regimes (quantum and 
classical). We then give an analytic treatment of quantum mechanical refiection in the 
thin wall regime, in which interactions with the plasma can be neglected as a particle 
propagates across the wall. We focus on leptons because of their much weaker coupling to 
the plasma. We argue that they are likely to be accurately described by this calculation, 
but quarks are not. The relative magnitude of the baryon asymmetry produced for different 
fermions depends on their relative Yukawa couplings {not their zero temperature masses), 
their transport properties and their interactions. We calculate the baryon asymmetry for 
various parameter ranges and conclude that asymmetries comparable with observations 
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1. Introduction 

Almost twenty years passed before it was realized that the conditions stated by 
Sakharov for the dynamical generation of the baryon asymmetry in the universe were 
actually satisfied in the standard Weinberg-Salam model [1] . The required baryon number 
violating processes arise from the chiral nature of the theory and the topology of the vac- 
uum, as encoded in the chiral SU (2) l anomaly [2] . C and CP violation are built into the 
standard model, and Sakharov's final ingredient, a departure from thermal equilibrium, is 
believed to be provided by the first-order nature of the electroweak phase transition. 

Despite the fact that the minimal standard model has all of Sakharov's ingredients, 
it has proven difficult to use it to produce the required asymmetry. The main obstacle to 
its success is the very small CP violation in the KM matrix. A second problem is ensuring 
that the baryon number violating processes are suppressed after the phase transition - this 
translates into a rather low upper bound on the Higgs mass M//— 35— 80 GeV, potentially in 
conflict with the experimental lower bound. The exact constraint is still rather uncertain, 
because it relies on the details of the finite temperature effective potential, which has 
only recently begun to be be seriously studied [3]; [4], [5], [6], [7]. A further problem for 
baryogenesis in in the minimal standard model is that for the observed top quark mass, 
unless the Higgs is relatively heavy, M// — 130 GeV, the zero temperature effective potential 
is unbounded from below, and this leads to an unacceptable instability [8] . 

One notable recent attempt to elaborate a minimal standard model mechanism is that 
of Farrar and Shaposhnikov [9] , who make use of subtle finite temperature effects on the 
reflection of quarks from the bubble walls formed at the phase transition. This calculation 
has been argued in [10] and [11] to greatly overestimate the baryon asymmetry, because it 
neglects the imaginary part of the quark self-energy. A suggestion as to how CP violation 
could be amplifled in the minimal standard model was made by Nasser and Turok [12]. 
The idea here is that CP could be spontaneously broken by the formation of a Z boson 
condensate on the bubble wall, and that the competition between macroscopic domains 
could enhance the difference in free energy between the different sign condensates. 

The most clearly viable mechanisms elaborated to date work with extensions of the 
standard model, the simplest such extensions being those where there are extra Higgs flelds. 
From a particle physics standpoint this requires little justiflcation and readily provides the 
extra source of CP violation which seems necessary. Following suggestions of Shaposhnikov 
[13], and McLerran [14], Turok and Zadrozny [15] showed how this could work in multi- 
Higgs extensions of the standard model. In these theories there is a CP odd Higgs fleld 
phase which changes in a deflnite manner in bubble walls as the Higgs flelds 'roll' from the 
unbroken to the broken symmetry vacuum. This phase couples to the anomaly through a 
triangle diagram producing a biasing of the anomalous processes [16], [17]. As the bubble 
walls propagate through the universe they leave behind a trail of baryons. 

Cohen, Kaplan and Nelson (CKN) proposed two potentially more efficient mechanisms 
in the same theories, 'spontaneous' baryogenesis [18], and the 'charge transport' mechanism 
[19]. The 'spontaneous' mechanism was an application of an earlier idea employing GUT 
baryon number violation [20] , the charge transport mechanism developed an idea involving 
a lepton number violating phase transition arranged to occur at the electroweak scale 
[21]. Dine et. al. [22] also gave a discussion of 'spontaneous' baryogenesis employing 
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electroweak baryon number violation in general terms involving theories with a higher 
energy scale. An interesting extension of these ideas is the work of Comelli et. al. who 
showed that finite temperature effects can cause the spontaneous violation of CP in the 
minimal supersymmetric standard model. Bubbles of both CP types nucleate, but any 
small explicit CP violation is greatly amplified because it comes into the exponent of the 
bubble nucleation rate [23] . 

CKN originally argued that 'spontaneous' baryogenesis would work well in multi-Higgs 
theories in the 'adiabatic' limit of slow and thick walls. A phase redefinition of the quark 
fields was used to identify a term in the Lagrangian which acts like a potential for fermionic 
hypercharge on the bubble wall where the relative phase of the Higgs fields is changing. 
CKN calculated the local thermal equilibrium established in this background, with cer- 
tain 'fast' perturbative processes allowed to equilibrate, and showed that it resulted in a 
non-zero baryon number violation rate. In [24] however we argued that inappropriate con- 
straints were imposed on global conserved charges in this calculation and that the thermal 
equilibrium which is in fact approached has zero baryon violation. A correct calculation of 
this effect involves finding the departure from this equilibrium induced by the motion of 
the wall, determining how the motion of the wall competes with transport processes which 
tend to restore it to equilibrium. In the limit of perfect particle transport we showed that 
the system approaches the thermal equilibrium with zero baryon production. Dine and 
Thomas [25] have also questioned whether it is correct to treat the system as described 
by a fermionic hypercharge potential, pointing out that this manifestly does not have the 
correct behaviour as the Higgs vev vanishes. 

The second mechanism suggested by CKN, the 'charge transport' mechanism, was 
argued to apply in the case of thin bubble walls. The reflection of fermions off the bubble 
walls in a CP violating way leads to the injection into the unbroken phase of particle 
asymmetries which then bias the baryon number violating processes in this region. In [21] 
this mechanism was applied to a heavy (4th generation) majoron neutrino with explicit 
CP violation in Yukawa couplings. Subsequently [21] CKN showed how in a similar way 
the effect of CP violation in the Higgs sector in a multi-Higgs extension of the standard 
model could be mediated into the unbroken phase by reflected top quarks, and sizeable 
baryon asymmetries produced. These calculations were supposed to apply only to thin 
walls because the 'quantum mechanical' reflection effect was taken to be suppressed when 
the wall is thicker (i) than the quark mean free path, because of scattering or (ii) than the 
Compton wavelength, because of the validity in this limit of the WKB approximation. 

As pointed out by CKN this second mechanism has a great advantage over previously 
suggested ones. It is 'non-local' in the sense that the CP violation and baryon number 
violation are separated from one another. The effect of the CP violation on the wall can 
be propagated away from the wall into the unbroken phase, where the baryon number 
violating processes are unsupressed. 

In this and an accompanying paper [26] we will consider these ideas in a framework 
which makes clear their relation to one another. Essentially there is just one physical 
mechanism: on the bubble walls a background held is turned on which affects both the 
dynamics and the interaction rates of the fermions. Our calculations will apply both to 
the case of a two doublet extension of the standard model and the case of a Z boson 
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condensate in the standard model. We will find that the property of non-locality can 
in fact be generalized to the case of thick walls and used to greatly enhance the baryon 
asymmetry in that limit. Corresponding to the reflection in the thin wall case there is 
a classical force on the wall which perturbs the particle densities on and in front of the 
wall, as well as a 'spontaneous' baryogenesis effect, albeit in a guise consistent with the 
criticisms of this mechanism in its original form. In the thin wall case we carry out new 
analytic calculations which take into account a problem in the original calculation pointed 
out in [24]. We develop a new analytic diffusion formalism in which all the distributions 
are determined dynamically. We concentrate on the case of leptons emphasizing that they 
are very efficiently transported into the unbroken phase and are likely to be best described 
by this calculation for typical wall thicknesses. 

The structure of the paper is as follows. In the next section we discuss B violation in 
the local thermal equilibrium situation which is taken to pertain around and on the wall 
during the phase transition and in particular we emphasize that total left-handed fermion 
number should be thought of as the driving force for B violation. In section 3 we discuss 
the general problem of the dynamics of particles in the background of a CP violating bubble 
wall. We discuss how the previously explored quantum mechanical reflection has in fact a 
very non-trivial classical limit. This suggests a quite different approach to the problem of 
modelling the response of the plasma to the CP violating background when scattering on 
the wall must be taken into account. The rest of the paper is then concerned with the thin 
wall regime and this other 'classical' calculation is deferred to the accompanying paper 
[26]. In section 4 we calculate the flux injected by reflection off the wall. In Section 5 we 
describe our calculations of the propagation of the injected flux into the unbroken phase, 
and give the derivation of a set of coupled equations for the particle species which describes 
the diffusion and decay of the injected asymmetries. In section 6 we solve these equations 
and calculate the baryon asymmetries produced in various regimes. In section 7 we discuss 
the issue of screening. While the diffusion formalism allows for a complete treatment of 
screening, this addition signiflcantly complicates the analysis, and so for the most part we 
ignore it. We attempt to justify this in Section 7, where we give a qualitative argument 
that including screening effects would alter our results at most by a factor of order unity. 
In section 8 we calculate the baryon to entropy ratio and discuss our results, comparing 
the cases of injected quark and lepton fluxes. In the concluding section we summarize very 
briefly and point out directions for further work. 

2. Baryon Number Violation 

The departure from thermal equilibrium which is required for dynamical baryogenesis 
occurs as the true vacuum bubbles sweep through the false vacuum. We shall in this 
and the accompanying paper describe this departure from equilibrium with space-time 
dependent chemical potentials for each particle species. Whilst the perturbations are not 
in general of this simple form, in many circumstances as we shall argue in detail below, 
thermalisation towards such a constrained local thermal equilibrium are rather efficient, 
and this allows for a greatly simplifled description of the baryogenesis process. 

The chemical potentials for the different particle species are fixed by the local values 
of the charges conserved by 'fast' processes which have time to equilibrate. Other 'slow' 
processes (such as baryon number violation) remain out of equilibrium over the relevant 
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timescales. In these circumstances the rate at which a charge X conserved by the 'fast' 
processes moves towards its equihbrium value is given by 



- AF 

X = -T^AX (1) 

where f is the equihbrium rate for the process which violates X, AF is the difference 
in free energy between the two states (in which the process is in and out of equilibrium 
respectively) and AX is the amount by which X changes in the process. This equation 
can be derived simply from detailed balance considerations. Applying this to the case of 
B violation one has simply 

B = -r.^^ (2) 

where Tg is the sphaleron rate (per unit volume per unit time)^ and //b is the chemical 
potential for B. CKN have phrased the problem in these terms and set out in each of their 
scenarios to show that //s is non-zero on or in front of the wall. 

It is useful to emphasize another form of this equation. Consider a process z^^nj ^ 
(in notation where, for example, 2ni +77-2 Sna has z/i = 2, = 1 and z/3 = —3) for which 
the rate is F. If the particles are in local thermal equilibrium with chemical potential //j 
then from (1) 

rii = - — {i^iT,jUjiij). (3) 

Using the fact that the sphaleron processes are tLtLbLTL ^ and tLbh^Li^L ^ (for a 
single fermion family with obvious notations) we can write find 

5 = - ^ (3/it^ + 3/i6^ + iir^ + 11^^ ) (4) 
or, if there are Np families, 

B = -NF^^iiSfiU + Sfil + fit, + fiij (5) 

where the sum is over the families. Here the chemical potentials stand for the difference 
between particle and anti-particle chemical potentials. Equation (5) tells us that what we 
must do in order to calculate B violation is follow how the densities of left-handed fermions 
and their anti-particles are perturbed. In their work on charge transport baryogenesis 
CKN have presented hypercharge or a charge X 'orthogonal' to it as the quantity driving 
baryogenesis. As discussed in [24] this emphasis is misplaced - it is clear from (5) that 
one can have non-zero hypercharge without any B violation, and B violation with zero 
hypercharge. In the approach we describe in these papers all charges will be determined 
dynamically and we will always use the formula (5) directly. 



't' We will use a bar over rates to denote rates per unit volume. 
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In the case that the fermions are massless there is a simple relation between the chem- 
ical potential and the perturbation to the number density: for particles Sni = 
and the opposite for antiparticles. In this case equation (5) becomes 

B = -QNf^{3Bl + Ll) (6) 

where Bl and L l are the total left-handed baryon and lepton number densities respectively. 
The factor of 3 is just a result of the definition of baryon number (i.e. ^ per quark) - the 
quantity in brackets is just the total left-handed fermion number. 

We can use (6) to comment on the role of strong sphaleron processes. Giudice and 
Shaposhnikov have pointed out [27] that if one treats these as 'fast' in the constraint 
calculations carried out by CKN (as they argue one should) one finds that the result is 
that there is no B violation. One must then go to the mass induced corrections to (6) 
to get a non-zero result. This result has a simple explanation. In all these calculations 
B and B — Ll are constrained to be zero. Also, since strong sphalerons couple right and 
left-handed baryons as tLtRbLbR.... it follows from (3) that the sum of the chemical 
potentials of left handed baryons is equal to that of right-handed baryons if this process 
is in thermal equilibrium. In the massless limit this leads to the relation Bl = Br which, 
together with the constraints B — — Ll , immediately gives B oc 3Bl + Ll — from (6). 
It is clear however that if one injects net left-handed lepton number into the system that 
this argument does not apply. We defer a full discussion of this issue to the appropriate 
point below. 

Now we turn to the analysis of how CP violation on the wall perturbs the left-handed 
fermion number on and in the vicinity of the wall. This then drives B violation according 
to equation (5). 

3. Particle Dynamics on a CP- violating Wall 

We are interested in the effects of CP violating condensates on the bubble wall. To 
date, the two simplest suggestions for this are the relative phase of two Higgs fields in 
a two-Higgs theory [15], [16], or the longitudinal Z boson [12]. Our analysis will include 
both of these possibilities. The former case is conceptually simpler - if there is explicit 
CP violation in the Higgs potential, the relative phase of the two Higgs fields is a CP 
odd field, which changes in a definite manner as the Higgs fields change from zero in the 
unbroken symmetry phase to their final nonzero values in the broken symmetry phase [15], 
[16]. The latter case relies instead on an instability which may occur if the top quark mass 
is sufficiently large [12]. 

Let us try first to define the relevant physical degrees of freedom. We can use the 
SU (2) symmetry to pick a gauge in which the background classical solution provided by 
the wall is (pi = -^(0, vie*^^)), (p2 = -^(0, V2e*^2)_ We assume that at all stages during 

the phase transition electromagnetic U{1) symmetry is unbroken. We could with a gauge 
choice set 9i or 02 to zero, but it is instructive not to do so, in order to clearly see how 
the (gauge invariant) CP-odd relative phase 9 = 9i — 02 emerges. It is clear that since the 
relative phase loses its meaning when the vevs vanish, all physical effects which depend on 
the relative phase of the two Higgs fields must vanish as either Higgs vev does. We shall 
see that this condition is met in a somewhat nontrivial manner. 
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The physical degrees of freedom in the broken phase are fluctuations in the two vevs, 
Svi and Sv2, in the relative phase 9, and the charged Higgs fluctuations (p'^. Diagonalising 
the Higgs kinetic terms we find 

where v'^ = + and = gl + g^- We ignore for now the kinetic terms for vi and 
V2-, the charged W bosons and the charged Higgs fluctuations t/?"*" (we shall return to the 
latter below). is the usual (gauge variant) expression in terms of W'^ and B^. Both 
9 and Z^^ are gauge invariant, and we take them as our definition of the CP violating 
condensate field and the Z field respectively. 

Assume for simplicity that only couples to the fermions via Yukawa terms. The 
phase of the Higgs field 9i can be removed by performing a (or hypercharge Y - the 
two are equivalent since Q = -\-Y is aiv unbroken exact symmetry) rotation on the 
fermions, at the cost of introducing a coupling d^9iT^^ i.e. a pure gauge potential for 
into the fermion kinetic terms. Let us see what effect this has. We may combine the 
d^9iT^ potential with the coupling to the field to find 

ipil^id^ - igAZi,^^)ip - mtpip (8) 

where gA = +5'/4 for up- type quarks and (left-handed) neutrino, and gA = —g/^ for 
down-type quarks and leptons, and 

g,z, = g,zO'-\^^a^e (9) 

where qa = \g- The vector contribution from the potential does not appear in (8) 
for the following reason. We treat the wall as planar, and assume it has reached a static 
configuration in the wall rest frame, with the background scalar fields being functions 
only of and the field ~ (0,0, 0, Z(^)) being pure gauge. We can then remove the 
vector term using the remaining unbroken anomaly-free vector symmetries: A hypercharge 
rotation leaves a vector piece which couples each fermion in proportion to B — L, which 
can itself be removed by the appropriate B — L rotation. Precisely the charge qa as given 
in (8) is the linear combination [^(Ts — Y) + j(B — L)]g. The remaining axial term cannot 
be gauged away, and (as we shall see) has a real physical effect. 

Let us briefiy mention how the calculation is done in unitary gauge. In this gauge, 
the massive gauge bosons are just the original gauge bosons i.e. the Goldstone modes that 
would have been 'eaten' by the gauge bosons are set zero. This condition is just that the 
gauge current corresponding to each broken symmetry (with generator and correspond- 
ing gauge field VF^), give by dC/dW^ be set zero at zero W^. Then the Lagrangian C is 
quadratic in W^, and there is no need for any 'eating' to occur. In our case, this condition 
is just iip*{l — crz)d^(p + h.c. = 0, or vfd9i = — V25^2-Then Z^^ = Z^, and it is clear that 
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the fermion rotation by 9iT^ induces a potential of exactly d^OiT^ into the fermion kinetic 
terms. After dropping the vector coupling, this agrees with (9). 

In describing the rotations required to bring the Lagrangian to the form of (8) we 
did not consider the Higgs fluctuations. As remarked by Dine and Thomas [25], a 
rotation such as the one we performed to remove the phases from the fermion mass terms 
introduces some di dependence into the Higgs-fermion interaction terms. To remove this 
one has to perform a rotation on the Higgs fields, which introduces some dependence on 
d^6 in the Higgs kinetic terms. In the limit that the vevs vanish this rotation is just a 

rotation of the Higgs fields (or —Y when we perform the subsequent hypercharge 
rotation to remove the vector term in the fermionic Lagrangian), and the induced term 
simply a coupling to the field which makes up the remaining piece of the pure gauge 
potential. The charges corresponding to the form in (8) are and —g/2 for the charged 
and neutral Higgs components respectively, (so that qa charges which are just proportional 
to a linear combination of T"^ — y and B — L charges are conserved in any interaction). 
Let us consider a little further the Higgs sector when the vevs are not zero. The 

rotation on the fermions then induces a 9i dependence in, for example, the charged 
Higgs-fermion interaction term ^^tifin- We now define the charged Higgs fiuctuations in 
unitary gauge. With the notation {pi = [ip'l ,vi/ \/2) (i.e. absorbing the phase into fi), 
and similarly for ip2, one finds the unitary gauge condition is 

Vt*9lj-'^l + <^2'*^M^2 = df,(p^*Vi + df,(ppV2 (10) 

Now we write the kinetic terms for the Higgs fluctuations: 

K=\d^<ft\^ + \d^<pt\^<x ^^^^ 

We shall see how a potential for emerges in the Higgs sector in a certain approximation, 
more restrictive than in the case of the fermions. If we drop derivatives of \vi\ and \v2\, 
then we find upon substituting (10) for the first term in (11) that it is quadratic in 
and dij,92- Let us make the approximation of ignoring such terms. If we now rotate the 
fermions by e*^^-^ to remove the phase from the fermion mass term, we must rotate ip^ by 
the same rotation to remove the phase from the Higgs-fermion interactions, and (p2 by the 
same rotation to remove the phase from the Higgs-Higgs couplings. Setting (p'^ = /ie~*^^ 
(P2 = /2e~*^S we find the second term in (11) becomes 

K oc \d^flv2 - d^f^vi + id^e^{flv2 - /2*^;i)P ^ \d^x + ^X^M (12) 

where x = /i V2 — /l'*^! is the physical charged Higgs excitation. Thus, in this approx- 
imation, one really can describe the effects of d^6i as a potential for charged Higgs 
excitations (this was originally argued to be true by Dine and Thomas [25]). 

One other assumption we made in deriving (8) was that only a single Higgs field 
coupled to the fermions. If both Higgs fields couple to the fermions we can carry out the 
same procedure to remove all the phases from the Yukawa terms involving one Higgs field. 
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but are left then with the phase 9 in the Yukawa terms involving the second Higgs field. 

There are alternative ways of writing this Lagrangian - we can choose where to put the 
relative phase but cannot get rid of it and make all the mass terms real. In what follows we 
will always consider cases where at most a single Higgs field couples to any given fermion, 
so that in studying the dynamics of a given fermion in the CP violating background we 
can always take its Lagrangian to be of the form (8) with a real mass term. 

So what is the cff'cct on the dynamics of fermions of turning on an axially coupled pure 
gauge field in the presence of a mass? And how are the distributions of the particles in 
the plasma affected by this background? We first consider the dynamics of a free fermion 
described by the Lagrangian (8). Because of translation invariance in the direction perpen- 
dicular to the momentum is a constant of the motion. With an appropriate Lorentz 
transformation, this may be set equal to zero and the problem becomes one dimensional. 
With the ansatz ij) ~ exp —iEt, the Dirac equation can be broken into two equations, one 
coupling the first and third components of •0, and a second coupling the other two: 

We are working here in the chiral representation to follow the convention of CKN in [21]. 
The energy E here equals ^/E^'^^^ in the original Lorentz frame. What (13) describes is 
the coupling of ingoing left-handed fermions to outgoing right-handed ones and vice versa. 
For anti-particles the same applies, except that the signs are switched i.e. left-handed 
particles and their anti-particles (which are right-handed) see opposite signs for the axial 
field. Since left-handed fermions and their anti-particles are affected oppositely we have in 
principle a way of generating a disturbance in left handed fermion number, which should 
source baryon number as discussed in section 2. 

This is precisely the effect identified by CKN in the 'charge transport' mechanism. 
Once one notes that the refiection from the wall is CP violating one must simply calculate 
the injected fluxes and determine how they bias baryon number violation. Emphasis has 
thus been placed on the calculation of quantum mechanical refiection and transmission 
asymmetries. We will discuss in detail in section 8 the criterion for the validity of these 
calculations in which the fermions are treated as free particles in their interaction with 
the wall. Roughly the criterion is that the thickness of the wall be much less than 
the mean free path Xmfp of the fermions. As the wall becomes thicker L„ ~ ^mfp strong 
scattering effects must be taken into account. Because CP violating reflection has been 
understood, for reasons which will be explained below, as a quantum mechanical effect it 
has been assumed that it is washed out in the limit L^j >> Xmfp- In the next section we 
discuss these reflection calculations, pointing out that there is a non-trivial WKB result 
which leads to quite a different conclusion and points the way to a new calculation of the 
thick wall case in which scattering of the fermions as they cross the wall is taken into 
account. 

3.1 Two Calculations of Reflection Coeflicients 

As discussed in section 2 the quantity of interest in to baryon production is the dif- 
ference in particle minus anti-particle distributions. We thus focus on calculating the 
difference in the reflection probabilities of particles of a given chirality and their antiparti- 



9 



cles, since this is the quantity that will enter the calculation (in section 4) of the currects 
injected into the unbroken phase. 

(i) Case 1 : Quantum Mechanical 

In Appendix A we derive an expression for the difference TZ in the reflection coefficients 
(i.e. probabilities) of particles and anti-particles. For left-handed particles L the result is 

n{p,) = Rl^r - R-L^R = \ / Im[m{z)m}] cos{2p,z) dz (14) 

where ruf is m(— oo), the mass of the fermion in the broken phase, and t — tanh^ where 
tanh2^ = |m/|/|p2|. Here Pz denotes the value of the momentum at infinity (in the 
unbroken phase). Henceforth we shall set the phase of the fermion mass m/ to be zero in 
the broken phase. Equation (14) is valid for any wall profile and gives the leading term in 
an expansion in L^nif, where is the thickness of the wall. Only when this expansion is 
valid is there a range of incident momenta which are both not totally reflected {pz > mf) 
and also non-WKB [pg < L~^). It is because of this that we refer to it as the 'quantum 
mechanical expansion'. 

To see the effect of the cosine term in the integral we evaluate (14) for an imaginary 
mass of a Gaussian form : Xm[m] = m/ = ^ exp — (mn-z)^- The result for the reflection 
coefficient is then: 

n = 4:t(l-t'')^ecpe~^^y , (15) 
rriH 

where Qcp is defined by J Im [m]dz — -^Qcp- (15) illustrates that the effect is quantum 
mechanical, as it is exponentially suppressed in the WKB limit, and in the limit \pz\Lyj < 1 
clearly involves a non-local coherence effect across the wall in which the particle and anti- 
particle pick up different phases as they propagate. 

This is the effect identified by CKN when they calculated the refiection coefficients 
numerically for a definite wall profile, and also in [28]. They found strong suppression of 
reflection for momenta larger than the inverse wall thickness and concluded correctly that 
they were seeing a phase coherence effect which would be strongly suppressed for all mo- 
menta when the effects of scatterings on the bubble wall are accounted for. Consideration 
of the WKB limit will now show that this analysis misses an important point. 

(ii) Case 2 : WKB Regime 

We now consider the limit in which the functions m and Z vary slowly in comparison 
to the Compton wavelength of the fermion. The WKB solutions to (13) then describe 
excitations which may be thought of as particles with definite momentum and position. 
They are described by the dispersion relation 



E = 



'pl + {^M+^T9AZf]'^^ S' = ±l (16) 



where refers to the spin of the particle in the frame where p± vanishes. For particles 
incident from the unbroken phase this is equivalent to chirality. This equation only makes 
sense for m > \gAZ\, and we assume this holds everywhere. If this condition is violated. 
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'positive energy' solutions can mix with 'negative energy' ones and one expects particle 
creation to occur (as in the Klein paradox). Wc shall not consider this possibility here. 
A particle incident from the unbroken phase with momentum p^(oo) = Poo is reflected 

if 

\Poo\ < m{z) T 9aZ{z) (17) 
holds for any z. From (17) we can read off the difference in the reflection coefficients: 

Rl^r - Ri^R = I 1' b°o I G [m/, max[m(z) + qaZ {z)]] ^^g^ 
~^ 1 otherwise 

(where we have assumed gAZ{z) is positive), ignoring barrier penetration effects. 

What this analysis reveals is that high (i.e. WKB) particles see the axial field 
on the wall as producing an extra potential superimposed on the real mass barrier. This 
extra potential has an opposite sign for particles and anti-particles. If the shape of the 
extra potential is monotonic like that shown in Figure 1 this potential does not produce a 
difference in reflection coefficients for the 'classical' particles and their anti-particles. The 
only interesting reflection will come at low 'quantum mechanical' momenta for which the 
different phases picked up by particles and antiparticles moving in this potential causes 
differential reflection. However if the potential looks more like that in Figure 2 (with a 
bump) there will be a difference in reflection coefficients of particles and anti-particles even 
in the WKB limit. Thus we see that the exponential suppression we observed in (15) was 
a result of the particular ansatz we took for the wall profile. Similar monotonic ansatzes 
have been assumed by CKN and other authors [28]. 

Besides revealing a sensitivity to the ansatz in the calculation of the reflection coeffi- 
cient, in particular that there can be non-trivial reflection of WKB particles, this discussion 
is of relevance to the thick wall case. The effect being treated was expected to be sup- 
pressed when scattering on the wall becomes important since it relies on a phase coherence 
of the particles across the wall. However the effect does not arise from such a phase coher- 
ence in the case of WKB momenta particles. As we have argued the picture appropriate to 
these momenta is of a (local) classical potential which is simply different for particles and 
anti-particles. What do we expect such a potential to give rise to? Consider the effect of 
turning on an electromagnetic potential in a small region of a plasma. Given a sufficient 
amount of time the potential will be screened - it will draw in charge - and settle down to a 
new thermal equilibrium in which there is a net overdensity of every particle in proportion 
to its electrical charge times the local value of the potential. If we consider the case in 
which the region moves slowly with velocity v the system will remain approximately in 
this thermal equilibrium but with deviations which vanish as f — > 0. We have just argued 
that the particles and anti-particles see an opposite potential and we thus expect that the 
deviations induced as the region - in this case the bubble wall - moves will reflect this 
difference. This should be true irrespective of whether the potential is monotonic or not. 
Just as in the case of an electromagnetic potential where the process of establishing the 
approximate equilibrium involves the electromagnetic force pulling in particles which then 
scatter in the potential and reach local equilibrium, there should be a force playing an 
analagous role. This force will be a CP violating force which will source perturbations in 
front of the wall much as the reflected flux does in the thin wall case. 
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What is the force a particle feels moving on the wall? To answer this we return to 
the dispersion relation (16). The velocity of the WKB particle with energy E is the group 
velocity 

_ p. 1 

where pz is specified by (16). A simple measure of the effect of the background is the 
corresponding acceleration, which may be calculated using pz = —dzE: 

dvz _ 1 {m^ {QAZm^y ^ 

IT - "2 ^ ^ :^vF^ 

with E and p±_ constants of the motion. The CP violating effect goes to zero as the mass 
squared does, consistent with the fact that Z may be gauged away in this limit. 

In what limit is this force associated with the background likely to be important? We 
turn to this question in the accompanying paper [26], where wc develop a new formalism 
to describe the effect. An accurate treatment of the 'nonlocal' classical effect becomes 
possible for thick, slow walls, where a fluid approximation may be used in order to take 
particle scattering on the wall into account. But in the remainder of this paper we turn to 
the consideration of the opposite 'thin wall' quantum mechanical limit, taking the fermions 
to propagate as free particles in the background provided by the wall. 

4. Injected Fluxes in the Thin Wall Regime 

We first make a few comments about the physical parameters which determine the 
applicability of the calculations we are undertaking. In these papers we assume that the 
phase transition is first order and proceeds by bubble nucleation. This is supported by 
detailed study of the perturbative effective potential at least for rather light Higgs {ruH < 
70GeV) [29] and, more recently, by non-perturbative studies [7]. The parameters which 
are crucial to an accurate determination of the baryon asymmetry are the wall thickness 
Lyj and the wall velocity Vyj. Estimates of the bubble wall thickness from perturbative 
calculations in the standard model [30] indicate Lyj 2m^"^ ^ and ~ 0.1 — 1. A 

recent detailed study by one of us (T.P) and G. Moore [31] in the minimal standard model 
indicates that for Higgs masses in the range 30 — 70 GeV or so, ~ 0.4 and ~ 25/T. 
If non-perturbative effects play a significant role this result could be significantly altered; a 
non-perturbative condensate makes the phase transition stronger with more super-cooling 
and hence larger Vyj. Likewise in two- Higgs theories larger or smaller Vyj^s may be possible 
in different regions of parameter space. 

The other parameters which are crucial to accurate calculations are those character- 
izing the particles in the hot plasma. One of our objectives in this work has been to carry 
out detailed calculations of diffusion constants and decay rates rather than give rough 
estimates of these parameters. This is particularly important because the emphasis of our 
work is on transport and the transport properties of quarks and leptons are quite different. 
An estimate of the length scale over which we can treat the fermion as free is specified by 
the fermion damping rate 7 in a hot plasma [32], [10], [11]. For the slow quarks and leptons 
{p < gT) these are 

T T 3 T 

^ 2asT ^ — 7/ a^T ^ — 7^ -a^ tan^ 9^T ^ — (21) 
0.0 oU Z lU 
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at lOOGeV, where 9yj is the Weinberg angle, and the subscripts indicate quarks, left- 
handed leptons and right-handed leptons. For hard thermal particles [p ~ T) even though 
the damping rates are logarithmically enhanced: 7 ~ f?^ ln(l/(7^), the numerical value does 
not differ by much: 7^ ~ AagT/S, ~ a^jT. Comparing these with L^, we see that 
the uncertainty in translates to an uncertainty about whether the fermions should 
be treated as free or interacting particles on the wall. It seems most likely with current 
understanding of the relevant parameters that leptons may be well described by a free 
particle treatment on the wall, but for quarks an interacting fluid approximation is needed. 
In this paper we shall consider the former case, in the companion paper [26] the latter. 

In the context of these remarks about the likely difference between the quarks and 
leptons we note an interesting point related to the finite temperature properties of a two 
Higgs doublet model, which we believe has been overlooked in the literature. The ratio of 
quark and lepton masses is not in general determined by their zero temperature ratio. For 
example the ratio of the top quark mass to the r-lepton mass at finite temperature in the 
models with the couplings described above is 

There is no reason why the ratio of the zero temperature masses should not arise from a 
tuning of the vevfs, in the potential rather than from a tuning of the Yukawa couplings. If 
there is such a tuning we do not expect it to be respected by finite temperature corrections 
to the potential. Rather we expect all the finite temperature vem to be of order T. Thus 
it is the ratio of the Yukawa couplings which is important and these can be large or small, 
consistently with present experimental data. 

We now turn to the calculation of the injected flux in the thin wall regime which 
we will use in the next section to calculate the departure from thermal equilibrium in 
the unbroken phase. In order to get an expression for the asymmetry injected into the 
unbroken phase, we need to include both reflection from and tranmission through the 
barrier. One can use CPT invariance to find the following relations between the reflection 
and tranmission coefficients [21] 

Rr^l = ^L^R = 1 - Tl^l = 1 - (23) 

where the transmission coefficients are for particles incident from the broken phase. In- 
tegrating the reflection and transmission coefficients against the incident flux, and using 
(23) to eliminate the transmission coefficients, we arrive at an expression for the fiiix of 
injected right handed particles 

where fll_ and (/:^^) are the free particle phase-space densities (in the wall frame) of left- 
moving left-handed (L) and right-moving right-handed (R) particles, respectively. When 
the wall is at rest the term — /j^ cancels so that there is no net reflection. Taking this 
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term to linear order in Vyj, and using the coefficients in (14), calculating again to leading 
order in and also to leading order in we find 



J/' = -JP = ^mjniHecP . (25) 

For the reflection coefficient we have taken 

nip,) = \ -jr^^CP mf<p,< rriH (26) 
L otherwise 

In imposing the upper cut-off we assume that there is a sharp WKB suppression like that 
in (15). The approximation 2t ~ tanh29 — ^ is good to leading order in since the 
integrals are dominated by the larger momenta because of phase space factors. 
In the WKB case, using (18), we obtain 

Jl^ = -J^ = ^m)mHQcpri (27) 

where we took gfA-^max = tuhQcp- The parameter r] parametrises our ignorance about 
the exact profile and penetrability of the barrier. For a thick (impenetrable) barrier with 
a prominent bump such that max[m(^) + qaZ] = nif + ^A-^maxj rj = 1. In general rj < 1. 
All masses in these formulae are the relevant finite temperature masses. 

What is the correct ansatz for the wall profile? Although monotonic ansatzes of the 
type previously assumed seem plausible they are no more than that. This point is well 
illustrated by the calculations in [33] in which it is shown that in the standard model the 
backreaction of the CP violating refiection can cause an instability to the formation of a 
Z condensate. Without a detailed calculation of these effects it remains unclear exactly 
what the wall profile is. 

5. Propagation of Injected Particle Asymmetries 

We now consider the problem of how the fiux calculated to be injected in the first 
two cases above propagates in the unbroken phase and drives B violation. The average 
velocity Vi relative to the wall of the injected flirx can be written as 

_ !p,<o4^^^^^~ ^^^'^^''^^^ (28) 



Using (14) and the same approximations as for the calculation of the fliix (to leading order 
in Vwi mf/niH and niH/T) we find 

- ' (29) 



' 4/n2 T 

For the WKB case, using (18) and (24), the result is 



Vi = 



(30) 



18C3 T 
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where Cs — 1-202 is a Riemann zeta function. The result in both cases is roughly the ratio 
of the momentum in the direction of the wall of the typical reflected particle to the average 
thermal momentum in that direction. Note that Vi is independent of the wall velocity. The 
particles in this flux propagate away from the wall until they scatter. If the wall is not 
moving too fast they will then thermalize and diffuse until they are overtaken by the wall. 
Any given particle can also decay into others through decay processes and these must be 
taken into account. 

For what range of wall velocity Vw do particles have time to thermalize before they 
are recaptured by the wall? Suppose that r is the mean time for a particle's velocity to be 
randomized. We take r to be both a thermalization time and the step time in an isotropic 
random walk which the particle executes once it thermalizes. Then we can estimate the 
mean distance a particle moves away in the direction of the wall motion in time t to be 

-y/^. Equating this to the distance v^jt the wall moves in the same time we see that the 

the ratio of the 'injection time' r to the 'diffusion time' t is approximately 3v^. We will 
restict ourselves to the case when this ratio is small i.e. 

«Vs = -j= (31) 

where Vs is the speed of sound in the relativistic plasma. We can then model the prob- 
lem analytically with equations describing the diffusion and decay of a particle density 
perturbation sourced by an injected flux localized in a small region around the wall. 

5.1 Calculation of Diffusion Constants 

We present our calculation of the relevant diffusion constants in Appendix C. We 
generalize the standard treatment of this problem given in [34]. We assume the local 
distribution function f{x,p) to be specified by a space dependent chemical potential plus 
a perturbation which we fix by balancing the contribution to the change in / due to the 
gradient in the chemical potential against the collision term in Boltzmann's Equation. 
For the dominant gauge boson exchange diagrams we then arrive at expressions for the 
diffusion constant D which defines the relation between the diffusion current and the 
gradient driving it: 

Jd = -DVn (32) 

where n is the number density. We will give our numerical results at the appropriate point 
below. 

5.2 Derivation of Propagation Equations 

We now use the continuity equation to derive the equations which describe the prop- 
agation into the unbroken phase. This gives 

^ = ^ + V • J. = - E ^-^(E ^^)- (33) 

A 3 

The decay term is derived using (3) and the sum over ^ is a sum over all the decay channels 
of particle i. Taking the leading (massless) term we can express this in terms of the particle 



15 



densities. In Appendix D we outline the calculation of the specific decay rates which we 
will need below. 

Now making use of (32) and (33) and treating the wall as providing a flux J/"^-' so 
that Jj = —DiVui + J/"'^, we flnd 

- a,n, - J2 ^ff^) = ^^J^r'- (34) 

A j 

The wall can be treated in this way since it conserves particle number, simply sending 
equal and opposite particle number in opposite directions. It can thus be described as 
producing a 'blip' in the flux, and taking it to be constant in the wall frame, we model it 
as 

jr\z) = ^a!S{z-v^t) (35) 

where Jf is the net reflected flux for each species z, given for example for right handed par- 
ticles by equation (24) . We deflne a parameter here called which deflnes the persistence 
length of the current in the vicinity of the wall, and we approximate the injected current 
with a delta function. This is reasonable if the injected current thermalizes in a time Tth 
short in comparison to the time a particle spends diffusing before being recaptured by the 
wall - which is just the criterion (31). 

Our ignorance about how the injected flux thermalizes is parametrized by ^j. The 
uncertainty in this parameter is unfortunately intrinsic to the analytic approximation we 
are using in taking the departure from thermal equilibrium to be modelled by a chemical 
potential, a form to which the injected fliix does not in general conform. 

Nevertheless an estimate of can be made as follows. For a diffusing particle with 
diffusion constant D the velocity randomization time r is ~ 6D. This can be obtained 
from the following simple consideration. The average distance a particle with a velocity 
randomization time r diffuses in a given direction in time t is (2^) = ^tr. On the other 
hand the solution to diffusion equation dtn = D'V'^n with a point like source: n{ f, t) ~ 
_|_g-r /4Dt specifies: {z'^) — (r^)/3 — 2Dt. A simple comparison of the two gives r ~ 6D. 
Using this decay time for a fiux injected at velocity Vi, we estimate ^ ~ QDiVi. ^ can also 
be taken to include other suppressions of the injected fiux which may be relevant. We will 
return to a discussion of this question in section 8. 

Finally we look for solutions which are stationary in the rest frame of the wall, i.e. of 
the form n{z, t) = n{z — Vyjt). The equations are then 

Diu'l + VyjTi'i - ^ ^a(^ kfrij) = ^iJi5'{z - Vyjt) ^gg^ 
A j 

where kf is the factor which results when we express the decay terms in terms of the particle 
densities and F^ = a. The problem now becomes that of determining the baryon 
density behind the wall (in the broken phase where we take the clcctroweak sphaleron 
processes to be turned off) when a flux of some species is injected in front of the wall. 

6. Solution of the Propagation Equations 
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To understand the solution of these equations let us first consider the case in which a 
chiral tau lepton flux is injected by the wall. The only decay processes we will consider are 
the Higgs mediated processes shown in Figure 3 and the electroweak sphaleron processes. 
In the unbroken phase there are no other processes which cause a chiral lepton to decay; 
in the broken phase there are other vev-suppressed processes but we will neglect these, 
taking them to be slow in a sense which will become more precise through our analysis. 
For the moment we will assume that the electroweak sphalerons are too slow to significantly 
alter the densities in front of the wall and thus that they too can be neglected, except for 
their role in making baryons. The Higgs processes bring about a coupling between the 
tau leptons and several other species ( Higgs particles, weak gauge bosons, quarks). For 
simplicity let us assume that these induced densities in species other than tau leptons 
can be neglected. We will later relax these assumptions and treat the full set of coupled 
equations. We then have 

DlL'[ + VyjL'^ - TLR{aLL - bLn) = ClJqS' ^ ^ 

(37) 

DrL'Ii + v^L'^ + TLR{aLL - bin) = -UJqS' 

where we keep a — ^ and 6 = 1 as variables for heuristic purposes. Left handed parti- 
cles come as isospin doublets, but only one particle participates in each interaction. This 
explains A;^ = a = 1/2. The problem of solving (37) reduces to the corresponding homo- 
geneous equations subject to the boundary conditions 

(38) 

DrLrI^ = -^rJq DrL'j^ + v^Lnlt = 

which are derived by integrating up the equations through z = 0, imposing the conditions 
that the and Lji are at most step- like discontinuous across the wall and that they are 
zero at -|-oo (in the unbroken phase). 

Substituting the ansatz e~^^^ into (37) we find that the exponents are the roots of the 
equation 

DlDri^^ - v^{Dl + Dr)ii^ + [vl - VLR{aDR + 6I?l)]/x2 + v^{a + h)VLRlJ. = 0. (39) 

This gives a constant solution and a cubic which has two real positive roots and one 
negative root. Requiring that the solutions be zero at -|-oo and finite at — oo we have eight 
constants to determine in the solutions: 

Ll = Lie-^'^ + 126-^"^^ z>0 

= Lse-^'-'^ + L4 z<0 

(40) 

LR = Rie-^''' + R2e-^''' z>0 
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where //s is negative. The conditions (38) fix only four of these constants. The 
additional condition required is found by adding the equations and integrating. This tells 
us that 

DlL'^ + DrL'j^ + VujLl + VyjLn (41) 

is a constant of motion. By fixing it everywhere we can determine all the constants. 

This conserved quantity (41) is simply the dynamical version of lepton number conser- 
vation which is built into the equations. More generally there is exactly one such relation 
for every quantity conserved by the decay processes we include in (36). B — L conservation, 
for example, will take the form 

D,,B'l + D,^B'^ + D^^T'^ - Di^L'^ - Di^L'^ + v^B - L) ^ (42) 

where is the diffusion constant for left-handed quarks etc. If all the diffusion constants 
are equal and the injected B — L is zero then B — L will be zero everywhere. However since 
in general the various species diffuse differently this will result in non-zero B — L locally, 
consistent with the constraint of global B — L conservation embodied in (42). 

If we can solve the equations (37) to find Ll we can then use (6) to find the induced 
baryon number. Taking the baryons to be thermalized they appear as a flux —v^B in the 
wall frame and therefore 

B = -v^d.B = -GNf^Ll = -T,Ll. (43) 

This equation is simply the limit of the appropriate diffusion/decay equation for baryon 
number which results when we use the fact that the diffusion constant of the quarks is 
very much less than the diffusion constants of the leptons, which source the equation. 
Integrating up Ll as in (40) we find the baryon number on the wall to be 

Thus we need simply to determine Li and L2 by solving the cubic equation for the roots 
and the boundary problem. 

We can solve (39) approximately in two different limits: 

Case 1. vl » TLRiaDR + bDL). 

In this limit the roots are 

Ii3 = -{a + b) ^i = — + + (45) 

Vw Dl Vw Dr Vw 

and, using the boundary conditions (38) and (41), we find 

Ll = -tt-Jq R2 = — -pp- Jo 

^ (da) 
, hVLRDR „ oTlrDl ^ ' 

L'2 — 9 til ti\ — ^ 

V V 
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From (44) we then have 



m--i^Jo)[—, (47) 

-t^ w w w 

We have taken ^l/D^ = Cr/Dr, assuming ~ Dj from our estimate above. The first 
factor in (47) is the amphtude of the diffusing flux, down hy ^l/Dl relative to the injected 
flux amplitude due to back diffusion across the wall. The expression inside the brackets is 
a conversion factor of this injected flux to baryons. 



;2 

Now the roots are 



Case 2. < « TLRiaDn + BOl) 



(49) 



1^1 = aDl + WL '''" = ^ ^^LR{aDl^+hD-^^) (48) 

and the solutions 

a ' aDn + bDL^DR aDR + bDL^DL ^' 

bPR gDr g L»L. 

aDii + bVL Dr aUR + 6L>i: Dl Dr 

to leading order in both v^/TlrDr and v^/TlrDl- (48) and (49) give 

B(0) = (^J„)££(^|^^. (50) 

a + b 

These two limits of these equations have a simple explanation. The mean time a 
particle with diffusion constant D spends in the unbroken phase before it is captured by 
the wall is ~ D/v'^. If, as in case 1, this is long in comparison with the decay time 
then the solution describes the two injected fluxes diffusing with a perturbation to each 
of order TlrD/v'^ which tells us how the two species are sourced by one another through 
the decay process. Likewise since B violation has been assumed to be a slow process and 
there is no injected B the conversion factor is for the same reason TgD/v'^. (47) therefore 
simply shows how B is sourced directly by the injected Ll and indirectly by the decay 
of the injected Lr into Ll. The precise coefficients can be read off from the decay term. 
We note that the signs of the contributions are opposite and therefore that there can be 
(for certain decay and diffusion constants) a critical wall velocity for which they will 
precisely cancel. 

The second case above is that in which the particles have time to decay before they 
are caught by the wall. In the ni solution which extends furthest into the unbroken phase 
the left and right-handed leptons are in abundances such that the decay process is approx- 
imately (to corrections of order v'^/TlrD ) in equilibrium i.e. Li = ^Ri- Working to 
leading order in Dl/Dr << 1 we have -Ri = ~^^'^o- The right-handed lepton density 
is simply equal to the effective injected flux and then putting the decay process in equi- 
librium determines the density of left-handed leptons. We can also read off fii from the 
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lepton number conservation constraint (41), since each exponential solution individually 
will satisfy this, by simply putting in the equilibrium constraint. The profile of left-handed 
leptons which act as a source in the baryon number violation equation (43) is 

--(-^"^0)6-^^^ (51) 

We note that this is different from what one would infer if one took the local densities to 
be determined by the equilibrium constraint and a lepton density taken to be carried by 
the right-handed leptons diffusing in front of the wall. This calculation gives the profile 

^ (|^</o)e-^^ (52) 



a^-h^Dn 

which, however, when integrated gives the same baryon number -B(O). 

We can in fact generalize the simple argument used to determine this sort of 'long tail' 
solution when there are n coupled equations, with one species which diffuses more efficiently 
than any other. Firstly we enforce the equilibrium constraints on all the appropriate decay 
processes. The equations then have n — A conservation laws of the form of (41) and (42), 
where A is the number of decay processes in equilibrium. We can write all but one of 
these as linear combinations excluding the species which diffuses in the long tail solution. 
It follows from the conservation laws that these linear combinations can be set equal to 
zero to leading order in the ratio of the diffusion constant of any other species to the 
one diffusing in the long tail. We now can determine all amplitudes in the solution in 
terms of the amplitude of the diffusing quantity. By substituting these back in the final 
conservation law involving this species we can then determine /ii for the solution. Finally 
we need to convince ourselves that the amplitude in this exponential solution really is equal 
(to leading order) to the effective injected flux for that species. To see this one writes down 
the boundary condition giving the jump in the density at the wall for the species which 
diffuses in the long tail e.g. 

Dr{Ri + R2 + R3 + ...) = ^Jo. (53) 

Using the same conservation law applied to all the other solutions individually one can see 
that the amplitudes of all of them (i?2,-R3 etc.) are next to leading order relative to the 
amplitude of the long tail provided the amplitude of the source for all other species is 
not larger than that for the diffusing species. 

Let us now summarize this procedure which we will use to extract the information we 
need to calculate the baryon asymmetry from the full set of decay equations : 

• 1. Identify how the injected fluxes of species with diffusion constants Di can source 
left-handed fermion number signiflcantly for the relevant wall velocity v^j. Either they do 
so directly or through some decay process with rate F which satisfles > 1. 

•2. Consider only the single species with the longest diffusion length Do which is 
injected (flux J^) and sources left-handed fermion number. The amplitudes we need to 
determine are those in a solution ni{z) = ni{0)e~^^ which describes the diffusion of this 
species and its sourcing of other species through decay processes {/i ~ v^/Dq). 
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•3. Identify all the A decay processes which satisfy > 1. Impose the condi- 

tion that the appropriate signed sum of the particle density amplitudes (assumed linearly 
proportional to the chemical potentials) corresponding to each process is equal to zero. 

•4. Identify the remaining n — A linear combinations of particle densities conserved by 
the A equilibrated processes in step 3, and write them so that the density of the sourcing 
species appears in only one linear combination. Set all but this latter linear combination 
to zero and so determine the amplitude of all species nj(0) in terms of the sourcing species 
amplitude, which is just jf^Jo (Co the appropriate persistence length). 

•5. To determine the root /i of the exponential decay tail corresponding to the solution 
use the remaining conserved linear combination in the form (42) with all diffusion constants 
but Do set equal to zero, and all the amplitudes in the coefficient of solved in terms of 
the sourcing species amplitude from step 4. 

•6. If the baryon number violating processes were amongst the equilibrated processes 
in step 3 the final baryon asymmetry is simply the amplitude of the baryon number in the 
solutions just determined. If the baryon number violating processes are not equilibrated, 
the solution for the total left-handed fermion number must be used as a source in B = 
-Ts{3Bl + Ll) which simply gives S(0) = -£^(35^ + Ll)(0). 

For our analysis of the case of injected leptons this will be adequate to determine the 
baryon number in all the different cases of interest. For top quarks the solution of these 
equations is complicated by the fact that the leading contribution from these 'long tail' 
solutions vanishes in various cases so that the steps following the second above are not 
applicable. We will then have to extend our analysis to find the contributions from other 
solutions which we will have to determine more carefully. However, as we will discuss below, 
we will limit ourselves in this case to a qualitative analysis of the resulting asymmetry. 

6.1 Baryon Asymmetry from an Injected Tau Lepton Flux 

We now consider in detail the baryon asymmetry produced in various regimes by an 
injected flux in tau leptons. First we note the numerical values of the parameters which we 
need. In Appendix C we calculate the diffusion constants taking into account the dominant 
t channel gauge boson exchange diagrams in Figure 5 

The results are 

= 8a2^(l + 0.8tan4^^)T ^ 
D-"- = 2%al, tan^ d^T ^ ^ (54) 
D-^ = SalT ^ I 

where 9w is the Weinberg angle [35] . Note that the right-handed leptons which are coupled 
only through hypercharge interactions have a very long diffusion length: Dji k, 380/T 
4Dl. 

In Appendix D we calculate the rates in a plasma for the perturbative decay processes 
mediated by Higgs bosons as in Figures 3. For these we find 

F, = 0.2a,yf Vlr = = 2Tr^ = 0.3a„y^ (55) 

where yt and i/r are the Yukawa couplings for the top quark and tau lepton respectively. 
The rates for the Higgs exchange process in Figure 4 is also calculated in Appendix D, 
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but we will not make use of it as it is subdominant for the cases of interest to us. The 
normalizations of these rates are defined by the convention in (36). Note that this means 
we do not incorporate the counting factor (for color and isospin) into the rate, just as in 
the case we have discussed where F lr was defined so that a factor of 1/2 appeared in front 
Ll because of the two isospin states. Ta in (36) is just the collision rate per unit volume, 
divided by a factor T^/12. 

For the anomalous processes we have Vi in (36) differing from unity and the convention 
we have adopted is an awkward one. Instead we define the rates as the coefficients in 

S = -r,(3Si + Lz,) A = -r,,A (56) 

where /S. = — B^. The second equation is derived by the same arguments outlined in 
section 2. These give 

Ts = QNFnsawT ^'^^IQ'^KsT Tss^Q^KssCi%T ^ Kss^- (57) 

where we have adopted the standard conventions in which KgOt^^T^ and l^sgCt^T^ are 
the number of topology changing processes per unit volume per unit time, and Np = 3 is 
the number of families. 

Numerical simulations give values of Hg and Kss in the range 0.1 — 1 [36] 

Case 1: vl > TrDR, T.Dr. 

This is simply case 1 as treated above with -B(O) as in (47). The Higgs processes do 
not have time to equilibrate the injected particles with any other species, and the baryons 
are produced directly from the injected left-handed Icptons and from the fraction of the 
injected right-handed leptons converted to left-handed leptons. 

Case 2: TrDn > vl > TsDr. 

Now the Ft- processes can equilibrate in the right-handed lepton diffusion tail. This 
creates a source for Higgs particles which couple in turn to quarks. We work for definiteness 
in a model in which one Higgs doublet (pi is coupled to the charge | quarks and the second 
doublet (^2 to the charge — | quarks and leptons. We take all the processes in Figures 3 and 
4 to be in equilibrium as well as the strong sphaleron. Other permutations and possible 
Higgs couplings can also be considered in a way which will be clear from this example. 
Following the general procedure outlined above we must therefore impose 
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where tl is the density of left-handed tau leptons and tau neutrinos (minus their anti- 
particles), 01 the density of Higgs particles of both charges in doublet 01, ql the density 
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of left-handed tops and bottoms, ul the density of left-handed quarks in the lighter two 
generations, etc. The numerical factors are the counting factors for color and isospin, 
and also a factor of 2 for fermions relative to bosons which enters in the conversion from 
chemical potential to number density. The gauge bosons cancel out because when we 
subtract the equilibrium constraints for the processes involving all the anti-particles from 
those for the processes involving the particles. The conservation relations which we use to 
solve for all densities in terms of tr are 

qL + tR + bR = 

TL + h-bR- \nL = 

1 (59) 
(pi-tR- -ul = 

nL + nR = 0. 

We find these by simply determining directly the linear combinations conserved by the 
processes which we have put in equilibrium. Together (58) and (59) give 3Bl + Ll = 
Ql + + tl = tl = ^^tr, and following the procedure outlined above we calculate that 
in the 'long tail' solution the profile of 3Bl + Ll is 

182 , , 300 Vw_ y 

:(|^Jo)e"^^ . (60) 



which integrates to give 



118^L>fl 



T}((\\ ^^'^^sDr , ^R . , . 
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We note that the correction obtained over a calculation neglecting all species except the 
tau leptons is negligible: the factor ^ replaces |, which follows from tl — 2tr and 
Ail = 3vyjD^^. 

Case 3: T^Dr, TsDr > vl > TsDl 

In this case the sphalerons have time to reach equilibrium in the Dr tail. We must 
now simply add the constraint 

qL + nL + TL + lL = (62) 

where II are the left-handed leptons in the two lighter generations. We also revise the 
conservations laws (59) by adding II with the appropriate coefficient to each of them so 
that they also invariant in electroweak sphaleron processes. We then proceed as before 
and calculate the equilibrium density of B at the wall (which is just the amplitude of the 
solution) to be 

Neglecting all species but the tau leptons we would solve 3Bl + tl — 0, Bl — Br and 
Tl = '2tr to get -B(O) = ^{-^Jo), which differs by a factor of two from what we have 
calculated including all species. 
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Case 4: T^Dl > vl > T^Dr. 

Here the tau lepton-Higgs processes are too slow to allow the right-handed particles 
act as a source for baryon number, but the left-handed leptons are diffusing long enough 
to allow the baryon number violating processes reach equilibrium. Since these couple the 
leptons to the quarks we must impose all the constraints in (58) above except the first one. 
We use the same conservation laws as in case 3, simply removing tl in the appropriate 
way from the second one. Solving for B{0) we find 

B(0) = -A.,(0) = (64) 

which again is a very minor alteration of the naive result B{0) = —^tl{0), which follows 
from imposing = Br and 3Bl + tl = 0. 

6.2 Baryon Asymmetry from an Injected Top Quark Flux 
We will not treat the case of an injected top quark flux in great detail. Unless the 
walls are very thin they are likely to be very poorly described by our analytic calculation 
because (i) we require << m.f.p, and (ii) we calculated all our expressions in the case 
that ruf < rriH-, which means for the top quark L^, ~ rn^ < m^^ ~ 2/T. We will discuss 
these criteria in further detail in section 8. 

In the accompanying paper [26] we develop a new formalism for the treatment of 
quarks when scattering is taken into acount. The case of top quarks in the thin wall 
regime has also been studied by CKN in their original exposition of the 'charge transport' 
mechanism in [21], and our main interest here has been to explore the interesting possibility 
that leptons might source this sort of baryogenesis. For the purposes of comparison with 
these previous calculations and with our calculations in [26] however it is of interest to 
turn to this case with the insight we have gained by studying the case of leptons. The 
analysis we sketch can be carried through rigorously and leads to the results which we give 
below. 

CKN calculated the baryon asymmetry by considering the injection of a hypercharge 
flux into the unbroken phase. The diffusion of this flux was modeled numerically and 
the resulting profile used as a source for a local thermal equilibrium calculation, in which 
the various conserved global charges were set to zero. In [24] we pointed out that these 
constraints are not strictly appropriate as the local densities of global charges may float, 
i.e. are not necessarily confined to their globally constrained values. In the case of leptons 
we have just treated this is very clear e.g. B — L becomes locally non-zero despite the 
fact that it is zero in the injected flux. The method we have developed above allows us to 
treat this dynamical determination of the global charges and it is interesting to ask what 
cff'cct this has in the case of quarks. A second correction to the calculation of CKN is 
that which is required because of strong splialerons. It was pointed out in [27] that CKN's 
local thermal equilibrium calculation gives a null result when these processes are put in 
equilibrium. The method we have developed allows us to see the effect of a flnite strong 
sphaleron rate and calcualte the correction it leads to. 

We will assume again that we have a model in which one doublet is coupled to the 
charge | quarks and the other to the charge — quarks and leptons. We will take the 
Yukawa couplings for the latter to be small so that the only Higgs coupling which is relevant 
over the time in which particles are diffusing is the top quark one. 



24 



The simplest case we can treat is when > TDq for aU decay processes. Then the 
injected left-handed baryon number drives B violation directly, and is caught before there 
is any significant coupling to any other species. We read off 

m = -^-^ipo). (65) 

where ^ is as before the persistence length of the injected current. 

However from (54), (55) and (57) we see that the equilibrium condition v"^ < TDq 
is actually satisfied for typical non-relativistic velocities ~ 0.1 for both the top/Higgs 

2 

processes (f^ < ^) and the strong sphalerons (f^ < ■^). We cannot proceed however 
precisely as we did in the analysis of an injected lepton flux for two reasons: 

• If we take the left and right-handed quarks to have the same diffusion constant the 
conservation law (in the absence of weak sphalerons) 

DqB' + v^B = (66) 

for baryon density B forces S = everywhere. In the 'long tail' solution in which we impose 
the equilibrium constraints on the strong sphalerons we have B^ = B^ and therefore, as 
explained in section 2, since Bl = Ll = 0, there is no force driving B violation. 

• Our analysis relied on the assumption that the sourced species which diffuses furthest 
is directly sourced by the injected current. This is not true now because the Higgs particles 
which are sourced indirectly by the decay processes have a diffusion constant L><^ ~ -Dl 
(for the left-handed lepton), much larger than that of the directly sourced quarks. 

To disentangle the implications of these two points we first neglect the Higgs processes 
completely. The first point is true irrespective of whether the Higgs processes are turned 
on or not, and when we have clarified its implications we will return to this second point. 
We consider then the equations for the left and right-handed baryons to which our diffusion 
equations reduce in the absence of the Higgs processes: 

Dq.B'l + v^B'l -^{Bl- Br) - 3T,Bl = ClJoS' 

^ P (67) 

Dq^B'^ + V^B'r + ^{Bl- Br) = UJoS' 

where and Dg^ are the diffusion constants for left- and right-handed quarks respec- 
tively. We consider these equations in two cases: 

Case 1: = Dq^^ = Dq. We add and subtract the two equations (67) to get two 
equations for B and A = Bl — Br 

sr 

DqB" + v^B' - -^(5 + A) = 

op ^ (68) 

DqA" + ^„A' - + A) - r,,(A) = 2^ Jo5'. 
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The Tg term can be neglected in the second equation, and the solution for A obtained 
is Jo exp — ^/T\^JD^z to leading order in v^,/ ^/V^Dq- This is then used as a source 
term in the first equation for B, and the resulting baryon asymmetry at the wall is 

We see that the effect of including the strong sphalerons is to produce a suppression by 
the factor v^i/'^.^V ssDq ~ Vyj/.Ty/K^g of the result in (65). 

The factor of | can be simply understood. It arises when we solve the second equation 
in (68). The boundary conditions impose a jump across the wall of^/D. When vl > TssDq 
the amplitude in front of the wall is ~ ^/-D and the amplitude behind suppressed. When 
< TssDq the amplitude of the solution behind is equal and opposite to that in front, 
and the roots equal, so that exactly half the injected fiiix appears in front of the wall, and 
half behind. The explanation for this behavior is simply that in the low velocity limit the 
decays play the role of dissipating the injected asymmetry so that the motion of the wall 
becomes unimportant and the solution approaches the symmetrical one. In the other limit 
the injected flux is dissipated by the wall catching up with the diffusing particles and the 
only stationary solution is reached by particles piling up in front of the wall so that the 
back diffusion current exactly cancels the injected current. 

Case 2: - Dg^ = SD 0. 

In this case the argument that there is no contribution from the long tail breaks down 
because B is not zero locally in front of the wall. To treat this case we consider the altered 
version of (68), in which there is now a source for B, look for the contribution from a 
solution with A = 0. The equation for B can be shown to be 

DgB" + v^B' - = -^^W (70) 

where we have taken ^ oc D for the two chiralities. The solution for the baryon density 
behind the wall is 



Dq l/12afT is the average of the two diffusion constants. 

If we take 5D to arise from the difference in left and right-handed quark diffusion 
because of SU (2) processes like those in Figure 5, we use our calculations for the case of 
leptons to estimate =- ^ {aw/ 0:3)"^ ~ 1/20. A difference in the diffusion constants also 

results when one incorporates the Higgs processes in Figure 3, for which the scattering 
rates differ by a factor of two for the quarks of different handedness. A naive estimate 
from our calculated rates (55) indicates that this would give dD of the same order. So 
for low wall velocities and/or large Kss (71) indicates that the dominant effect can still 
come from the diffusion tail. The suppression due to the strong sphalerons in (69) is cut 
off when one takes into account the slightly different transport properties of the left- and 
right-handed baryons. 
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Finally we return to the second difference we noted between the case of quarks and 
that of the leptons - that the particle species which diffuses furthest are the Higgs particles 
which we have assumed not to be directly sourced by reflection from the wall. Including 
the Higgs processes should not alter significantly the solutions we have just considered. 
As in the case of the leptons these processes will simply cause a redistribution of the 
injected asymmetries between the species coupled through these processes, and this will 
lead simply to minor numerical corrections. The important point is that these processes 
do not drive the quantity driving baryon number to zero as the strong sphalerons do. In 
neglecting the Higgs processes, however, we did overlook the new 'long-tail' solutions for 
the Higgs particles which are sourced indirectly through the decay processes. Can these 
particles which diffuse much further than the quarks themselves in turn source significant 
baryogenesis? By examining the boundary conditions in the way we did in the lepton 
case (choosing this time to write down a conservation law involving ^ and a combination 
of quarks which is sourced) one can show that the amplitude of this solution is down by 
approximately Dg/Dfj,, relative to the directly injected diffusion tail. When integrated 
therefore, we expect, the two solutions should give the same order contribution. However 
the further suppression due to strong sphalerons which we have just discussed will be 
greater (by ~ ^jDq/D(f)) for the longer tail since the baryons are diffusing longer, and 
therefore we anticipate that no significant additional baryogenesis will occur because of 
the efficient diffusion of Higgs particles. This qualitative analysis can be carried through 
rigorously just as in the lepton case by fully solving the diffusion/decay equations with the 
appropriate boundary conditions. 

7. Screening 

In the treatment of the problem we have presented we have so far entirely neglected 
the effect of hypercharge screening. The original mechanism of this sort proposed by CKN 
described the problem in terms of a hypercharge ffirx generated by the reflection of top 
quarks off the wall. It was pointed out by Khlebnikov [37] that this treatment overlooked 
the fact that any such hypercharge density would in fact be very efficiently screened by 
the plasma. Since it is hypercharge that was described as driving the production of the 
asymmetry it appeared that this will cause a very significant attenuation of the effect. 
This criticism applies to the mechanism as we have treated it as well If one computes the 
hypercharge density profile in front wall in any of the stationary solutions to our diffusion 
equations the result is non-zero. CKN in a later paper [38] argued that the effect of 
screening could be accounted for by doing their constraint calculation with injected y in a 
different basis of charges 'orthogonal' to Y . One of the charges named X which is a linear 
combination of hypercharge and baryon number is identified as the appropriate injected 
charge which drives the B violation. Here we will attempt to clarify this question of the 
role of hypercharge screening, justifying the approximation we have made in neglecting it 
and outlining how it can be incorporated in these calculations. 

Consider first a system in a constrained thermal equilibrium, with chemical potentials 
jjii for species i. Suppose now we apply a local potential for hypercharge (/'r(a:) in some 
region. The local thermal equilibrium with the same chemical potentials is 

fi(p,x) = — T? — ; — , , . ^ where = \/ -\- m?. (72) 
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This makes the colhsion integral on the right hand side of the Boltzmann equation zero if 
the original chemical potentials fii did. This is true because shifting the chemical potential 
by an amount yiCpyix) does not alter the rates of processes which conserve hypercharge. 
This distribution function describes a solution in which the potential has the effect of 
drawing in an amount of each species proportional to its hypercharge, but this is the only 
change admitted to the local density of any species. It is this local equilibrium which CKN 
calculated in [38], fixing the potential (pvix) by requiring it to screen the hypercharge 
density in (72) to zero. The same fii as in the unscreened case are obtained by their 
prescription of working with 'orthogonal' charges. These are simply the linear combinations 
of charges which are unchanged by screening i.e. charges Qa = S^q'^ni in which (prix) 
cancels out because Tr{QY) = Hikiq^yi = 0, where ki is a statistical factor which is 
one for fermions and two for bosons in the massless approximation. Imposing the same 
constraints given by a calculated injected fiux on these charges amounts to solving precisely 
the same set of linear equations as in the unscreened case, except for the one setting the 
hypercharge to zero. Quite simply one is taking the previous solution and superimposing 
on it the extra screening densities which leave rates unchanged. In particular one finds 
that the sphaleron rate is unchanged and hence the resulting asymmetry. 

Is this a good approximation to the real process of screening? It describes correctly a 
static situation where some real physical constraint — an injected decaying fiux or a force 
like that we describe in [26] — imposes the constraints on the /Xj in some fixed region. 
Then this static equilibrium (72) satisfies the Boltzmann equation. However the cases we 
are considering are not manifestly like this. The region in which there are perturbations 
induced (described by locally varying chemical potentials jjii) is moving and the screening 
is dynamical. The screening currents respond to the injected currents and the time it takes 
for them to screen the charges is important. In particular the static approximation breaks 
down because the response of different species to the driving hypercharge field is different 
because they have very different transport properties. To see this suppose a single species 
is injected by refiection from the wall. This injected fiux carries positive hypercharge one 
way and compensating hypercharge the other way, which induces hypercharge currents in 
the other species across the wall. If all these other species are very "stiff" in comparison to 
the original species i.e. have much smaller diffusion constants, the screening will impede 
the diffusion we have described. However, if there are species with comparable transport 
properties (as there are in the electroweak plasma), the injected fiux can be screened by 
these species and the effect should be simply that the injected fiux drags along screening 
charges with it, with its behaviour little altered from that we have described. This is the 
picture which we will now briefiy support with some simple quantitative arguments. 

The way to include this dynamics is simply to incorporate screening in our calculations. 
We do this by adding the number current induced by the hypercharge field to the diffusion 
current using Ohm's Law so that we have 

J = Ja + Jg = -DVn + aE (73) 

where E is the hypercharge field induced by the local densities and a is the conductivity. 
Note that this definition of conductivity differs from the usual one by a hypercharge factor 
J/, because it relates the number current rather than the charge current to the field. The 
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conductivity may be related to the diffusion constant using the derivation of the diffusion 
constant given in Appendix C (see also [24]). The one difference is that the space varying 
chemical potential is replaced by the hypercharge potential. One arrives at a relation 
between the induced number current and the electric field which is exactly the same except 
as that between the diffusion current and the gradient in the chemical potential, with 



a - ^yD (74) 

D 



Incorporating this extra term in the derivation of the diffusion/decay equations (36) simply 
results, after using Gauss' law, in an extra term —A-KCiY in each equation. So we have a 
set of equations exactly like those we analyzed except for this extra coupling between all 
the species, which greatly complicates the analysis we presented. 

In the limit in which all diffusion constants are equal, a general analysis is possible. 
Taking linear combinations we can rewrite the N equations (for N species) as N — 1 equa- 
tions in which the screening term cancels out, and one equation for the total hypercharge 
Y in which the decay terms cancel out (because total hypercharge is conserved). This 
corresponds precisely to going to the 'orthogonal' basis described by CKN. In this case 
then the only alteration to our equations is to add the screening term —ATiTnyiaiY (where 
the sum is over species) to the hypercharge propagation equation. This term is formally 
just like another decay process term and we can follow the sort of analysis given in section 
6. The "equilibrium" condition (now for F = 0) becomes (f^/47rSj(TiD) << 1. Thus for 

Vy, < ^yDT (where f = T^iyf ) the screening will be effective in the 'long-tail' diffusion 
solution, and in this regime we approach the static limit discussed by CKN. The new set 
of constraints on our 'long-tail' solution will be satisfied by superimposing a density of 
each species in proportion to its hypercharge on our previous solution. Note also that any 

'short-tail' solution in which Y will be restricted to a distance {■\J^yT)~^ of the wall, 
which is simply the inverse Debye mass in this model. 

The diffusion constants of different particles in the electroweak plasma are however, as 
we have seen, very different, varying over two orders of magnitude. The current induced in 
response to a hypercharge field is therefore not proportional simply to the hypercharge of 
the responding species. The question which interests us is which particle species screen a 
given injected species and whether this can have any very significant effects on the solutions 
we analysed. In particular we would like to see the effect of the very different transport 
properties of the quarks and leptons in the screening of injected flirxes of either. We will 
not carry out a general analysis of the equations but limit ourselves to a simplified set 
of equations to illustrate the intuitive argument given above. The most important point 
is that the "stiffness" of quarks will not prevent the diffusion of leptons, because there 
are other leptons around which will screen the injected flux. Rather than drawing along 
quantities of each species in proportion to their hypercharge, the leptons will be screened 
by other leptons, not by quarks. 
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So consider the following simple model: 



DaA" + v^A' - ATTGAiVAA + vbB + ycC) = saS' 

DbB" + v^B' - AiraBivAA + vbB + ycC) = (75) 

DcC" + v^C - ATTacivAA + vbB + ycC) = 0. 

The species A is sourced, while the species B and C screen. We have not included a decay 
term as this should not be important because all such decays are hypercharge conserving. 
Let us now study the roots of these equations and compare them to the unscreened analogue 
in which A simply diffuses as described by its own diffusion constant (i.e. in the solution 
with root Vw/Da)- 

To obtain an equation for the root of the screened solution we combine (75) to get 
two equations without the screening term which can be integrated once to give two first 
order equations. Requiring e~^^^ to be a solution and imposing Y — gives 

DADBDc{J2y^)>l 

(76) 



— V. 

+ 



I 

vn, [y\DA{DB + Dc) + ylDB{Dc + Da) + ylDc{DA + Db)\\d 
vl[DAy\ + DByl + Dcyl] =0 



Solving this equation in the limit Dc « Da, Db we find a 'short' root involving Dc and 
a 'long' root 

, _ DAy\ + Duy% , for i^^ = i^^ = 

""'-'""DADBiyl + yD^X^:^^ iorDA»DB ^''^ 

We see that if there is one species C which has a very small diffusion constant, this does 
not shorten the previous solution provided there is a second species B available to screen 
which has comparable (or better) difi'usion properties. 

For baryogenesis we will typically want to evaluate the integral under the diffusion 
tail J A. In this case, to leading order in Dc/D, species C does not screen in the long tail. 
Using this and the fact that the main contribution comes from the longer tail, one obtains 

2 

f A = ^ y^+y^ ■ "^^^ suppression with respect to the unscreened case for which J A = ^ 
is due to the screening with species B. The second limit in (77) shows that the diffusion 
solution for A does not survive if both species have much shorter diffusion lengths — the 
solution is forced to follow the behaviour of the least 'stiff' of the screening species. 

We now solve (76) in the limit Da ^ Db « Dc- We assume Da = Db = D which 
considerably simplifies potentially cumbersome algebra but still models the relevant effects 
of screening. This case is more complicated because there are two comparable roots 

2 
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A careful analysis using the appropriate boundary conditions and imposing the condition 
of zero total hypcrcharge with the ansatz I = Is exp[~\sz] + 1^ exp[— Aj,^] for z > and 
1 = for z < (where I = A, B, C) permits us to evaluate the amplitudes. We find that 
in the shorter tail (A5) only species B screens; C does not screen at all (i.e. its amplitude 
is suppressed by D/Dc)- In the longer tail (Al) something rather unexpected happens: 
species B anti-screens and C screens such that yBB/ycC = —y%/{y\ + y\)- To answer 
how screening affects baryogenesis we again as above study J A. Of course, if species B 
and C carry left handed fermion number, the final baryon asymmetry should include terms 
proportional to J B and J C, but let us ignore this complication. Recall that if one species 
(B) screens with Da = -Ds = D we have J A = ^{yB/UA + Vb)- When one adds a 
second species (C) to screen with Dc » D one finds that the same result J A for the 
contribution of the shorter tail (in which C = 0). But the total result including the longer 
tail contribution is just J A = sa/v^ so that the unscreened result is recovered. 

In the electroweak plasma there are many species to screen, in particular there are 
lighter families which have almost identical diffusion properties to the third family fermions. 
For example, in the case of baryogenesis sourced by the right-handed lepton diffusion de- 
scribed in section 6, we expect that the screening will be provided by all three families of 
leptons (mainly by the antiparticles of the right-handed leptons). Still however it will be 
only the tr leptons which are converted to the left-handed leptons which source baryo- 
genesis and we thus estimate in the spirit of the simple model above that the net result 
would be an attenuation by a factor of |. The second case studied above tells us that 
something rather nice occurs in the top quark mediated baryogenesis. The screening with 
other quarks alone has the tendency to decrease somewhat the baryon asymmetry. How- 
ever the presence of mobile leptons and Higgs particles works in the direction of recovering 
the unscreened result. We conclude that the corrections due to screening of the results we 
calculated in section 6 should be numerical changes of order unity . 



We now calculate the baryon asymmetry in its standard form before comparing the 
various cases we have considered. Assuming the asymmetry at the wall to be frozen in as 
it passes behind the wall i.e. that the weak sphaleron is switched off just behind the wall, 
we need simply divide B{0) by the entropy density s = ^^g^T^ {g^ ^ 100 the number of 
relativistic degrees of freedom) to get the standard baryon number to entropy ratio. Using 
the flux in (25) we find, for the most interesting parameter ranges. 



8. Baryon to Entropy Ratio and Validity of Calculations 



TrDR >vl> TsDl 
TrDR, TsDr > vl 
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(79) 



^ A recent study of this question [39] has developed the present analysis further quan- 
titively and reaches the same conclusion. 
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where the first three cases are for leptons and the other two for quarks. In each case ^ and 
D apply to the appropriate fermion (left-handed lepton in first case, right-handed lepton in 
following two cases, quarks in final two cases). ©cP was defined by J Im [m\dz = ^^Qcp- 

As remarked earlier the mass m/ is the tree- level finite temperature mass, which we 
should take not to be the zero temperature mass but yf(j){T) ^ yfT where yf is the Yukawa 
coupling of the fermion. In particular for leptons in a two doublet model this parameter 
is not constrained to its standard model value by phenomenology [40]. 

The way we have written the result in (79) breaks it up into two pieces — an 'injected' 
piece on the left of the bracket and a 'conversion' factor on the right. It is clear that the 
leptons do much better in terms of conversion because of their much greater diffusion 
constants and the absence of the strong sphaleron suppression. The extra factor of three 
which the quarks gain is a color factor which really belongs on the left, because all three 
colors reflect off the barrier. The factor ended up in the conversion factor because we 
defined the injected flux for quarks in terms of baryon number. 

In section 3 we discussed how Z could be non-zero in either a two doublet model or 
in the standard model if a Z condensate is formed on the wall as described in [12]. In the 
former case it is appropriate to take Ocp to be simply the integrated phase change across 
a single bubble wall as this will be determined to be the same on every bubble wall by the 
effective potential, and we can take Ocp ~ 1 consistently with phenomenology. 

In the Z condensate case ©cp will contain some suppression (possibly many orders of 
magnitude) which will depend on the details of how one sign of the spontaneously formed 
condensate comes to dominate over the other as the transition is completed. This will be 
discussed in [33]. 

We note that for < TgD the asymmetry goes as 1/vw This is cut off if the weak 
anomalous processes have time to reach equilibrium. For velocities < (rg-D;^)^/^ ~ 
0.75y^q;u, cot^ Ow {cf. (54) and (57)) the asymmetry goes as Vw Note that for the right- 
handed leptons with their very long diffusion tail this means Vyj < O.ly^ so that the third 
expression above can be appropriate for modest velocities. For quarks on the other hand 
Vw < (FsDg)^/^ ^ 1.2y^a^/Q;s ~ requires velocities smaller than those in the 

calculated range Vw ~ 0.1 — 1. Recall that ~ 0.1 — 1. 

Before using (79) to calculate some numerical values, and in order to make a detailed 
comparison between the cases of an injected quark flirx and an injected lepton flux, we 
return to the conditions for the validity of (79). In writing it down wc have used the injected 
flux which is calculated from a reflection calculation which involved various assumptions. 
Firstly, we solved the Dirac equation for a free fermion by expanding it perturbatively in 
mfLyj. This expansion is valid provided 

• Condition 1. < mj^ (perturbative expansion) 

This condition can be avoided by solving the Dirac equation exactly for a given wall 
profile (Higgs mass profile) and then treating the imaginary part of the mass as a pertur- 
bation (see [28], [9]). For L^^, ^ ^/^i discussed in section 3, we expect a strong WKB 
suppression (for typical monotonic wall ansatzes, to which we restrict ourselves). 

Secondly, we treated the fermion as free on the wall. A typical particle in the injected 
asymmetry current has a momentum pg ^ Evz ~ mu ~ 2/Lu,, with typical energy 
E ~ 2T. The distance the particle then advances before it scatters once, given by VzlJ^, 
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should be larger than the wall thickness Lyj . This gives 
• Condition 2. Ly, < {l/T^f)^/^ (free fermion) 

One might also wonder whether a more stringent condition exists - namely whether 
collisions occur on a length scale shorter than the de Broglie wavelength of the particles. For 
particles of momenta ~ m/, the de Broglie wavelength is substantially larger than m^^, 
and a correct quantum treatment must include collisions. However, as noted above, the 
reflected asymmetric flux is dominated by over-barrier reflection of particles with momenta 
Pz ~ and for these Condition 2 suffices. 

We note that this condition is a stronger constraint than the naive ^ 7/ ^ which 
we used in section 4^. If condition 2 is not satisfied one can estimate the resultant phase 
space suppression. 

Particles with a large incident angle on the wall (p± Pz) are nonrelativistic so we 
can use pz ~ mfVz ~ "^jL^. The condition that these particles not scatter leads then to a 
weaker version of condition 2: Ly^ < (2/m/7/)^/^. A simple argument based on counting 
particles with large incident angle indicates that in this case an additional suppression of 
order l/{LwT)'^ to (79) will result. For the case of WKB reflection from a 'bump' barrier 
this condition is further relaxed because the reflection is dominated by thermal particles 
with momenta Pz ~ m/ . 

Consider now whether and how the quarks and the leptons meet these conditions. 
First consider the top quark. Since its mass is so large mt ~ T, the condition 1 requires 
very thin walls < 1/T, which is to be compared with the perturbative value L^, ~ 20/T. 
Next condition 2 reads: Lyj < 2/T. (The weaker version does not help since nit ~ T.) 
Hence unless the wall is very thin we expect additional strong suppressions for the top 
quark which have not been included in (79). 

For the leptons the conditions are much more plausibly satisfied. Take for example 
the right handed lepton whose long diffusion tail ends up in most cases dominating baryon 
production. The above conditions become: 1) Lyj < mj^; 2) Ly, < 9/T. This means that 
for a standard model Yukawa coupling (y^ = 0.01) and a wall Ly, ~ 10/T all conditions 

1/2 

are met. The weaker version of condition 2 gives: L^ < 9/mr so that for thicker walls 
we expect some phase space suppression. Note that we can allow y-r to vary quite alot 
(consistently with the phenomenology of two Higgs doublet models) without violating the 
conditions. 

In summary the differences between injected quark and lepton fluxes are: (i) leptons 
gain in 'conversion' in front of the wall because of their better transport properties and the 
absence of a strong sphaleron suppression which affects the quarks, (ii) for very thin walls 
~ 1/T the ratio of the injected asymmetries in quarks and leptons is equal to the ratio 
of their Yukawa couplings squared, (iii) for thicker walls ~ 10/T there are suppressions 
for the quarks which are absent for the leptons because of their much smaller mean free 
path; the treatment we have presented breaks down in this case, and new techniques for 
computing quantum mechanical reflection in the presence of scattering are required. 

Two other comments on (79) should be noted. The 'persistence length' ^ we estimated 
with a naive argument to be 6Dvi, where Vi is the injected velocity. The uncertainty in 



't' We are grateful to Larry McLerran for a discussion of this point. 
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^ reflects the limitation of our calculation. To calculate this parameter more precisely 
involves going beyond the diffusion approximation to determine exactly how the reflected 
asymmetry sources the diffusion equation. This involves modelling how the very specific 
momentum space distribution of the injected flux thermalizes as it moves away from the 
wall. The result is clearly ansatz dependent, as was manifest in our estimate for the 
quantum mechanical and WKB cases (for which the appropraite Vi were quite different). 
For a recent study of this question and a discussion of the limitations of the diffusion 
approximation see [41]. Finally we recall that our diffusion approximation relied on taking 
Vw < Vs ^ l/\/3: the speed of sound in the plasma. For highly relativistic velocities the 
particles do not have time to thermalize before they are caught by the wall, and clearly the 
asymmetry will be considerably attenuated. For > Vg the diffusion tail does not exist 
because particles cannot diffuse faster than Vg. In this case we expect local mechanisms to 
be the dominant sources of baryogcncsis. 

There is one aspect of the free particle approximation which we have not discussed 
- the neglect of the contribution of the thermal gauge boson excitations to the fermion 
self energies [42], [43]. Weldon gave a Dirac equation incorporating this and one can 
argue that this is the relevant equation for fermionic excitations if one is interested (as 
we are) in processes whose time scale is long compared to the response time of the gauge 
plasma ~ ^/gT, but short in comparison to the scales on which particles scatter, which 
is typically given by ~ l/g'^T or 1 / g"^ / g'^)T depending on particles' energy. The 
question of interest is then how the analysis of the reflection problem will be modifled. We 
present in Appendix B a set of manipulations of this flnite temperature Dirac equation 
which lead to a real space Dirac equation, the starting point for an analysis in terms of 
relection coefficients. We write the equation in a form in which one can see simply how the 
modification of the fermion self energies at finite temperature enters. We find that the mass 
matrix which relates left and right handed particles becomes off-diagonal as a consequence 
of thermal effects. These off-diagonal elements are of order q;^(LT)^ and hence small for 
a sufficiently thin wall. For realistic wall thicknesses however the corrections may be large 
and further analysis is required. 

Finally we turn to the numerical evaluation of (79) in a few cases, to illustrate 
that our final result gives asymmetries typically compatible with the observed value. 
For the standard model Yukawa couplings ruf/T Ki yf the prefactor in (79) reads 
g~^^v^{yf /L^T)'^QcPi we take ^ = 6Dvi, Vi as given in (29), m/ = yjT and 2mjj^ = 
Lyu. In the case of the r lepton the first formula applies even for rather slow walls: 

2 

Vyj > 2.1yr , \/kJ20. In this case we find ^ ^ -— /r\-,a — a^,,6cp- If we take 
y^ = 0.01 (standard model value), with L-wT ~ 20, g^f — 100 and Vyj ^ 0.1, we get 
ub/s ~ —0.6 X lO^^^KgOcp, marginally compatible with the nucleosynthesis bound 
(4 — 7) • lO"-^-*^, provided CP violation is large ©cp ~ 1- 

The result in other regimes can be read off simply. In view of our comments above one 
of the most interesting regimes seems to be that of intermediate Yukawa couplings y,- ~ 0.1, 
which is described by the second case above. The answer is larger by 3 x 10^ than the case 
of standard model Yukawa lepton coupling. In this case we expect the refiection calculation 
to be reliable, as the conditions discussed above still hold. An extra factor 3 was obtained 
from the conversion of the injected right-handed asymmetry to baryons. 
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9. Conclusion 

In this paper we have developed a new analytic formalism to treat the diffusion and 
decay of a chiral flux injected by reflection from a CP violating bubble wall. Our formalism 
allows us to fix dynamically all the injected quantum numbers rather than constrain them 
locally as was done in the previous calculation of this mechanism in the case of quarks. We 
have emphasized the case of leptons, in particular the interesting role that could be played 
by the right-handed leptons when the lepton Yukawa coupling is larger than its standard 
model value. We have also briefly examined the case of quark reflection, determining the 
suppression which results from strong anomalous processes. 

There are many interesting issues which remain outstanding some of which we have 
discussed briefly (See for example [44]). The injection of the current could be modeled in a 
more sophisticated way [41]. Our calculations apply in the limit when the wall velocity is 
not relativistic, but we anticipate they should be accurate until close to the speed of sound 
in the plasma. It would be interesting to understand the highly relativistic case in more 
detail too. The 'conversion' will clearly be very much less efficient but some of this may 
be compensated for by an enhancement of the reflected asymmetry as the moving wall will 
be Lorentz contracted. 

Appendix A. Quantum Mechanical Reflection Coeflicients 

We start from equation derived in the text 

from the Lagrangian (8), in the frame in which the wall is at rest (Z = Z{z) and m = m{z)) 
and in which the transverse momentum has been boosted away. E is the energy of 
the incident particle in this frame, related to the energy E in the unboosted frame by 
E — a/ E"^ — p^. We will consider particles incident on the wall (0 > 2; > ^o) from the 
unbroken phase on the right {z > 0). At z > zq we take the mass to have its (real) broken 
phase value mo. 

We first perform the global rotation 

iiaaSA Zdz 

C± ^ e ^± (81) 

so that the redefined ^± obey 

idz^± = Q±e± = ±1 (82) 

where 

i2igA I Zdz 

and m-t = me "^^o . This is precisely the form in which the problem is analyzed 
by CKN in [21], an analysis which we follow until (86) below. The solution of (80) can be 
written 

e±(z) = Pe"'/>±(^)'^^^±(^') (84) 
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where P indicates path ordering. The upper component of ^± is the left-moving eigenstate 
in the unbroken phase, so to describe the reflection of an incident left-handed (L) particle 
we put ^+(0) = (l) and ^+{zo) = iD~^ (g). D is the matrix which diagonalizes Q± = Qo 



in the broken phase, DQoD~^ = y & — m^a^. f and t are the reflection and tranmission 
amplitudes. We will analyze only this case, extracting the reflection coefficient for incident 
L (anti-particle of L) and incident right-handed particle R by the substitution m — > m*. 
Using these definitions and the explicit form of the solution (84) we have 

(J) ^ (f ) (86) 

where Q.{z) = Pe~^^o Q+{z)dz _ j^^^^-j _ (^*^' '^^^^^ ^i complex 

numbers to be determined. D can be obtained as ( * 1 where s = sinh6', sinh2^ = 
, and c = cosh^, cosh 2^ = , We then have 

sp -\- ca* sp -\- CO* 



We now define 



for which 



or 



U{^)=[^* a*)="'^^''" ^^^^ 



mle'^'^^ 



id,U= ( " "^+" ]U (88) 



za^a = m+e^'^'b* id,b = m+e'^'^^a* (89) 
which can be integrated recursively to give 



"Z . fZl 

a = l + 



f m+(^i)e2^^^i / 'm;(^2)e-'^^^^+.. 
JO Jo 

h=-i( I m+{zi)e^'^'^ + [ m+e^^^^i / ' m^e-^'-^^^ / ' m+e^^^^^ + 
\Jo -/o -/o Jo 



(90) 



Now using a = exp —iEzoa{zQ) and /3 = exp —iEzQb{zo), expanding the expressions (86) 
perturbatively in the integral dzm^{z) ex.p2iEz we find, to leading order. 



(^t + i £ m+e-2^^(---o) + it" m*e+2^^(---o)^ (91) 
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where t = tanh^ (tanh2^ = Squaring to find the reflection probabihty 



Rl^r = = - it{l - t^) 



20 



(izm+e-2»^(^-^o) 



= + t{l - t^) {2 j ne{m+)sm2iE{z- zq) -2 J Im{m+) cos2iE{z - zq)^ 

(92) 

Thus, using the replacement m m* for anti-particles, we obtain the leading term for the 
difference in reflection coefficients incident with momentum pz = ^^(-l-oo) = E: 



n{pz) = Rl^r - Rl- 



yR 



4t(l - t' 
|"^o| 



Xm[m{z)mQ] cos{2pzz) dz. 



(93) 



This formula should be treated with caution as at next order in the perturbative expansion 
there is a t-independent term 



2 I dzi / dz2'Im[m*_^_{zi)m+{z2)]sin2pz{zi — Z2). 
Jo Jo 



(94) 



Clearly the approximation (93) is only valid provided t is not so small that (93) is smaller 
than (94). In the cases we actually apply (93) this is only a marginal approximation. 

Appendix B. Finite temperature plasma effects 

The one loop finite temperature corrected dispersion relation for fermions in a hot 
gauge plasma was firstly discussed by Klimov in Ref. [42] . Here we follow Ref. [43] which 
gives a Dirac equation in momentum space which may be written 



[(1 + a)f+ 6^ - m] * = 



(95) 



where is the four-momentum of a particle and u^^ is the four-velocity of the plasma. 



In the plasma frame 



aL,R = 



MIr 

p2 



1-^ln 
2p 



(1,0), while in the wall frame = (7, 
cle. F( 

u;+p 



and m is the 



tree-level mass of the particle. For to <T , p < T the coefficients a and b read: 

,2 X 1 



LO — p 



>L,R 



MIr 



P 



+ 

. P 



p' 



- 1 



-In^^ 
2 LO — p 



(96) 



where Ml and Mji are the one-loop finite temperature corrected left- and right-handed 

1/2 

fermion masses, and uj = P^u^, p = \uj'^ — P^P^^ [43]. Since we are mainly interested 
in studying thermal particles with uj,p > T we have evaluated the integrals {A4) - {A5) in 
[43] and found that the relevant change to a and b is: ln(w +p) ^ l^C^^: where C ~ 1 — 2. 

The thermal masses for the left-handed (qL), right-handed up (uR) and right handed 
down (dR) quarks [43] are 



3 
8' 



= T^'^OigT^ + ^irayjT^ + — ttq;^ tan^ OwT^ 



1 

72 



M^ji = -TTasT"^ + -h -TTCKu, tan^ 9wT 
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(97) 



M^jj = -TTUgT^ + — TTCKu, tan^ OwT 
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while for the leptons 



= + ^««,T2 + tan^ OwT'' 

8 ^ 8 (98) 

Mf^ = + + -a^ tan^ 6'^^^^ 



For both quarks and leptons 



Ml-Ml^^ ^a^T^ (99) 



We incorporate the effect of a classical axial field in (95) by a canonical replacement: 

r^r-9A^6 (100) 

since the Z field couples axially. Wc wish to answer the question of how much the results 
for the reflection calculation deviate from the free particle case studied in the main text and 
in Appendix A. Our analysis will be far from exhaustive. We will content ourselves with 
showing that the corrected Dirac equation can be written with some simple assumptions in 
a form which is manifestly that of the free particle case, with corrections which are small. 

We consider as in the text a one-dimensional wall lying in the xy plane. In order 
to study the reflection problem we need to return to a real space Dirac equation. We 
do so by replacing —idz and making the following assumption: we assume a and 

b in (95) to be well approximated as scalar functions of u by taking p = p{uj) in (96). 
This is true in the leading order of a derivative expansion in the background fields i.e. in 
the WKB limit. Here we are however interested in the quantum-mechanical momenta Pg 
typically of precisely the inverse length scale of the variation in the background. Since 

p ~ [P]_ + P^~\ [P]_ — d1~\ this assumption should be valid provided p » 

which will be true for the thermal particles we wish to describe. As long as we restrict 
ourselves to particles incident at small glancing angles (which dominate phase space) our 
real 2;— space equation should apply. 

We start by re- writing (95) (with Z field included) for the two component right-handed 
"ifji and left-handed ^'^ spinors in the chiral representation (in which ^' ~ [^fR, ^l]) 

[idg - qaZ] = CT3 \{E + cru^) - cro • ^1 ^'jt + 02,mR{z)^L 
\iQz + 9aZ\ *l = [-(E + clw") - cl(T • *l - (T3mL(2;)*it 

where we used 75\I/i?,L = ^^r,l- As in the free case we are again working in the wall 
frame in which we have transformed away the momentum parallel to the wall by a 
Lorentz transformation 

P^=li.{P^-v^E) = ^ (102) 

such that 

^± = ^, 7± = -, E=Je^-PI (103) 
h, E ^ 
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and 



^ -> 

Ul_ = —^_i_V_i_U 



(104) 



Recall that = (0, 0, 0, Z{z)) and tt'^ = (k°, 0, 0, us) {u^ — 7^;, us — —'Jw'^w) are the field 
and the velocity of plasma in the wall frame. As a last step to obtaining (101) we have 
divided (95) by (1 + aL^ji) so that 



m{z) 
1 + ap L 



cr,l 



>R,L 



1 + Ofl,L 



(105) 



Note that (101) is more complicated than the corresponding free particle one since the 
chirality eigenstates \I/i?,L = (J) and '^r^l = (i) are no longer the eigenstates of the 
Dirac operator in the massless limit. Indeed the term 



-cr^losct ■ u 



(106) 



mixes the state (q) with (°). To find the true momentum eigenstates one diagonalizes 
the the Dirac operator: 



I^R = CTS 



{E + CrvP) - CrG ■ u 

-{E + cl-iP) - clB ■ u 



(107) 



This induces a nontrivial transformation on uir and rriL and they become non-diagonal 
matrices, i.e. a pure right-handed (left-handed) travelling wave couples via plasma effects 
to both components of the left-handed (right-handed) travelling wave (in the presence of 
a non-zero tree level mass). 

The momentum eigenvalues for the right-handed particles are 

= -CRUZ ± \r (108) 

and similarly for the left-handed fermions, with index R ^ L. It is trivial to re-write these 
relations in the plasma frame (in which u^^ — (1, 0, 0, 0)) 

= ± [{E + CRf -Pl]~' , P. = ± [{E + CLf - Pi] ^ (109) 
Note that the high momentum limit (P >> M) of these are 

^2 ^ p2 ^ 2M| , = p2 + 2M| (110) 

Roughly speaking thermal particles in plasma (P ~ T) behave as free particles with mass 
squared equal 2M^ ^ plus the tree level mass squared, with the important difference that 
plasma mass does not mix chiralities 



Ar = -CRU3 ± 



(E + CRU^) 



2 -2 
CrU± 
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Having found the eigenvalues (108) for Vji^l, it is now straightforward to construct 
the matrix that diagonahzes them 



7^ 



R 



E + cru^ 



(111) 



(112) 



where f^"* are the eigenvectors which correspond to 

Mr'' = 2{E + CRVP){E + CRVP + Xr) 

and similarly for the left-handed particles. We can now use TZr^l to diagonalize T>l,r'- 
W^VrUr = diag{-CRU^ + Xr, -cru° - Xr) = dR 

(113) 

HI^VlHl = diag(-CLW° - Al, -clu^ + Xl) = dL 

so that the left moving particles (negative eigenvalue) correspond in this new basis to ( ° ) 
for the right-handed particles and (q) for the left-handed particles and vice versa for the 
right movers (with positive eigenvalues). 
Now (101) can be recast as 

[id^ - qaZ] iiR = dR^R + mR{z)lCQ^L 

(114) 

[id^ + qaZ] = dL^L - mL{z))CQ ^^r 
where we have defined 7?.^ "^r = "^r. /Cq is an approximately diagonal mass matrix 

/Co = T^R^ crz T^L = Ko(T3 - rjasa ■ u± (115) 

Indeed one can check that 

«o = l + o(Ac)^ rj^-^^(l + o{c)) (116) 

2{E + c) 

In order to make a comparison with the free particle case we rotate the spinors in (114) 
one more time as follows: 

^R^^R^e'I'-'^'^" e"^^'-^R, ^L^^L = e~'f^^^'^' e"^^'^L (117) 
so that (114) reduce to 

id.^R = mR{z)e^'I'-^'''K:^L , id.^L = -mL{z)e-^' f '^^"'JC'^r (118) 
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where 

yv._e yv.oe -e l^r/e-2^Az(«i + ^2) -g-^^A^ ^ ^^^^> 

Where = {E-\-cY — c^u\. The matrix /C is diagonal for free particles (when T = 0) and 
equals 0-3. The deviation is due to thermal effects and it is of order 'qu±_ ~ Q!u,(TL)^. This 
is small for a sufficiently thin wall, such as we concentrate on here, but 0(1) or greater 
for thicker realistic walls. We leave a fuller analysis of the effects of these corrections to 
future work. 

Appendix C: Diffusion Constant Calculations 

To calculate the diffusion constant we generalize the method outlined in [45]. In the 
companion paper we present an alternative derivation based on a local thermal equilibrium 
ansatz for the distribution function, which yields a very similar result. 

We focus on the elastic scattering, t-channel vector boson exchange diagrams, which 
are expected to dominate the scattering process, and play the major role in limiting the 
diffusion of particles (the divergence in the propagator as t ~ M^, where M is the mass of 
the exchanged vector boson, leads to a logarithmic enhancement of the scattering rate). 

The diffusion constant D measures the particle number current induced by a gradient 
in the number density: 

= -DVn (120) 

where Jd is the diffusion (number density) current, n number density. Wc calculate D 
by considering a steady state situation in which some external source and sink cause the 
particle number density to vary slowly with x. To a first approximation, the distribution 
function is just /o = l/[exp(/3£'+^) + 1] (where /3 = 1/T is the inverse temperature, E = 
the energy, and ^ the chemical potential), and ^ varies slowly with x. This distribution 
function however has no current, so we compute the deviation 5f in the distribution 
function f — fo + Sf, hy solving the Boltzmann equation, and then compute the current 
from 

J = J <r'p5f I (121) 

The perturbation Sf is of order t/1 where r is the mean free path and I the length scale 
over which ^ varies: the calculation is accurate in the regime where r/l << 1. 
We determine 6f from the static Boltzmann equation, which reduces to 

— * 

dtf=^- V^f = -C{f) , (122) 

Expanding the collision integral C(/) in terms of 5f and expressing the L.H.S. of (122) in 
terms of V^n, we find the relation (120) from which wc read off the constant D. 

The collision integral for two in-coming {p^, k^} and two out-going particles {p'^, 

J{p',k,k'} (123) 

v[f] =fMi - fp'){i - fk') - fp'fk'ii - /p)(i - fk) 
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where S {^Pi) ensures the energy-momentum conservation, = J ^ pf^P^i and A4 
stands for the appropriate scattering amphtude. For particles exchanging B vector bosons 
(corresponding to the U{1)y symmetry), W bosons {SU{2)l) and g gluons {SU{3)c), the 
amplitudes are 

9 

\Mw\' ^ ^Qg L, _ (124) 



\M,\'^32gl + 



(t - M|)2 

where Y is the hypercharge of the particle we are following (in the convention Q = T^+Y), 
gi, g2 and gs are the U{1)y, SU{2)l and SU{3)c gauge couplings, with gi = g2tan9w, 
s = {p + kf = 2p-k, t = {p-p'f = ~2p-p', u = {p-k'f = -2p-k' and sin^^vK = 0.23, Ow 
is the Weinberg angle. Note that in calculating (124) we have set masses of all particles 
at the outer legs to zero. This is justified since we are interested in leading order {g^) 
contribution to the scattering amplitudes. The t-channel processes posess a logarithmic 
divergence when integrated over the momenta if the exchange particle is massless. One 
can cure this divergence by using the thermally corrected gluon propagator. Since in the 
approximation when the scattering particles are massless and on shell, t < 0, the gluon is 
Debye screened. To the leading-log accurate one can use for the propagator gr^,^/ [t — Mg), 
where Mg{T) is the Debye mass. In this approximation the dominant contribution to the 
scattering integral comes from the particles with the exchange momenta of order Mg. Note 
that so far we have not addressed the question of the difference between the longitudinal 
and transverse bosons. We suspect that the contribution from the transverse bosons, which 
are not Debye screened at the one-loop level, is larger, but since we do not know what 
is the approprate screening mass, we will treat all exchange bosons equally, as if they 
were longitudinal. The (one-loop) Debye masses of the (longitudinal) gauge bosons arc: 
Ml = ^a^ta.n^OwT'^ ~ 0.04T2, M^, = ^a^T^ ^ 0.69T2, = Syra.T^ ^ S.GT^ 
[43], where we wrote — 5'2/47r, as — 5'3/47r. Note that the gluon thermal mass is rather 
large: this means that the leading logarithm approximation calculation (which relies on 
an expansion in M/T) may be not very accurate. 

For the local equilibrium distribution functions, P[/o] = in (123) so that C(/) = 0, 
a consequence of the energy conservation p^ -{- k^ — p'^ + k'^ . Now we make our key 
approximation — we assume that there is little energy transferred in a collision, so that 
p ^ p\ which implies k k' . For light exchange bosons this is quite a good approximation. 
Nevertheless even for exchange of rather heavy gluons, this is true on average, since a 
particle is just as likely to gain as to lose energy in a collision. With this, we find to first 
order in 5/, V\f \ iSfp—6fp>)[f^{l — f^)), where we ignore the fluctuations in the second 
fluid 5fk, 5fk'. 

If locally particle concentration varies in x direction, then the L.H.S. of (122) becomes 
^da;^{x)f, where /' = df/d{pp^ + C). This motivates the Ansatz 



Sfp = gip,x)^ (125) 



42 



so that Sff = g{p') ^ ~ g{p) cos a, where cos a = p- p' /pp' ~ p ■ p' jp^ ■ 
The Boltzmann equation (122), (123) can be now recast as 

Jiv'.k.k'} 



2{pO) 



3 



Finally using (120), (121), and (126), with n ^ no - we find 

12 f ^3 fp 



The integration required for Tt is simple when one notices that the integrand of the k', 
and p' integrals is Lorentz invariant and can therefore be performed in the center of mass 
(CM) frame, in which p + k = = p' + k',p^ = k^j p'^ = k'^ . We compute to leading 
order in M/T, obtaining for example for the gluon exchange case 

r,.±Z!ln(iPlfi) (128) 

The diffusion constant follows immediately from (127): 

Ds' -^«^tan4eH.F^Tln(^ 



_ 1 80 /32T2~ 



/32T2\ 



keeping the leading order logarithms. (We evaluated the integrals containing logarithms 
using J x'^f'j, ln(l/a;) fa ln(l/a;„) / a;"/^, where fa n is the value which maximizes x^f'^. 
This is quite a good approximation because ln(l/a;) varies slowly where x^f'^ is large.) 
These considerations lead to the following diffusion constants for the right- and left-handed 
leptons and quarks 

D-^ =D7} fa — 
Ir b 

=D-i + 1 n-i ^ — (130) 

II W ' ^ Ir -^qq ^ ' 



D-,Ir =YLrK + 'l.rD^' + D-^ fa I 

where for the coupling constants we have used ct^o tan^ Qw — 1/103, ct^^ = 1/30, ag = 1/7, 
and for the logarithms involving Mb, Mw and Mq, the values 6.1, 3.8, 2.2 respectively, 
ez, = 1 and ei? = take care of the fact that the right handed quarks do not couple to the 
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W boson. To get a feeling how good is the expansion in (M/T)"^ we have also evaluated 
the next (mass independent) term and got a correction to the leading logarithm: 3/2, 3/2, 
9/4 for the S, W and gluon exchange respectively, which indicates that, except perhaps 
for the gluon case case, this expansion is a rather good one. 

Appendix D: Calculation of Decay Rates 

In this Appendix we outline the calculation of the decay rates due to fermion/Higgs 
and Higgs exchange processes shown in Figures 3 and 4. 

The rate per particle per unit time can be defined in the real time formalism as the 
integrated scattering amplitude squared in a thermal bath of particles: 



12 r 

J{P,P' 



fpfkil + fp'){l - fk')^\p + k-p'- k')\M\' (131) 

k,k'} 



Factor 12/T^ = fi/noT {p is the chemical potential and no number density for a single 
species) converts the rate per unit volume and time into the rate per unit time (see (32)). 
In the above J p = J c! p/2Ep\ fp and fk are the distribution functions for the m-going 
fermion and boson; (1 + fp') and (1 — fk') are the spin blocking for the oui-going particles. 
We assume thermal distributions for particles; if not thermal, the appropriate distribution 
functions should be used. 

First we focus on the calculation of the fermion/Higgs rate. The scattering amplitudes 
squared for an m-coming left-handed, right-handed lepton and a quark are: 



I A/f l"^ — A 



(132) 



where the constants are for a S-boson (on external leg) Ab = IGirawtatn'^ Owuf , for a 
VF-boson Aw = '^Anauiyf, and for a gluon Aq = ^TvasUq^ {6w = Weinberg angle, y's 
= Yukawa couplings, — 1/30, aw'ta.n'^Ow = 1/103). For fermion masses we take the 
one-loop thermal values: m^^ = 0.13T, m^^ = 0.33T, niq — 0.88T. To evaluate the integral 
in (131) the expression we need to integrate is 

12 r 
r I fpfkl 



T3 



^ (133) 

4. , , , ,x . -ts 



! ^"^ip + k-p' -k')A 

Jv'k' 



p',k' (t-m^)^ 

By ignoring the spin blocking we overestimate the rate by a factor 1 — 4. (The spin blocking 
has a value in [1/4, 1] for any momenta p, k, and it is approximately one (~ (1 — 2/e^)) 
for thermal momenta ~ 3T.) I can be easily evaluated in the CM frame (see Appendix 
C) 

A [ , 2p-k. 1 1 , , 

h(l + ^) - 1 + TTHlJ (134) 
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The integrals over p and k we approximate using: 



Vm{n) = J 



da;a;^e^^ In a; 



ml 



m+l 



n 



I 



1 1 ^ , ' 
1 + - + .. + C-lnn 

2 m 



1 



dxx^fx In — f=i In 



1 



a; 



(135) 



(C = 0.577 is the Euler constant, Xm+i ~ m + 1) and obtain the rate per unit time: 



647r3 



In 



(136) 



This estimate together with (132) gives the following rates for right-handed, left-handed 
leptons and quarks 



LR 



Ks 2 
2C3 



1 \ 

1 H — tan^ 9w ] In 

6 I m 



+ - tan ©w In 



4x2T^ 



m. 



T fti 0.28a^yfT 



37r2 



"syq In 2 



m 



(137) 

In order to obtain the approximate expressions for F's in (137) we have used tan^ Qiy = 
0.29, and the standard one-loop formula for thermal masses of fermions = g'^CnT^/S 
{Cr is the Casimir constant for representation R). This gives m/^^ = 0.13T, m/^ = 0.33T, 
ruq = 0.88T, so that the corresponding logarithms are ln(8T^/mf ) = 6.3, ln(8T^/m^ ) = 



4.3, ln(8T7m: 



2.3 {X2 ^ 2). 



The third family rates are fastest (r lepton and top quark) because of their large 
Yukawa couplings. 

In order to evaluate the Higgs exchange diagram (Fig. 4) we need the scattering 
amplitudes: 



M =A 



where A takes account of coupling constants and counting. Some care is required to count 
the diagrams for quarks and leptons. If one assumes that all fermions couple to the same 
Higgs (just like in the standard model) for the third family one gets 



At, 



--Ab, 



- ^'^yhh + 2(yt^ + yl)yl + other families 
—12y1y1 + ^y^yl + other families 
A^^ =l2ylyl + 4ylyl + other families 
A^^ =2A^^ = 12{yl + yl)yl + other families 



(139) 
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We now integrate (135) for the Higgs exchange process in the CM frame. The only- 
difference in the scattering amphtude is (— s) — > t: 



Stt 



m 



p ■ k 



In 1 + 



2p-k 



+ 



1 + 



2p-k 



(140) 



We can then approximate the integrals over p and k by expanding in [m/TY . The leading 
order term gives 

AT , , 

r = 141) 

512 - Stt ^ ' 

Taking account of counting (139) we finally obtain the rates for third family: 





=r6. = 


T 


1287r _ 
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2 2 1 

VtVb + 




1287r 
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2 2 , 

VtVh + 




1287r 




=2r., -- 


T 
1287r 



1 



ytyt + ^iyt + yDy^ 



.2„,2 



.2„.2 



+ other families 
+ other families 
+ other families 



(142) 



[iyt + yh)yl] + other families 



These rates are model dependent. Indeed if one takes perhaps by particle physics a more 
motivated case in which one Higgs couples to the up type of quarks and the other to the 
down type and leptons, the terms j/f j/^ and yly^ in (142) drop, altering the rate. 
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Figure Captions 

Figure 1: Classical potential seen by a particle crossing a bubble wall. The solid line shows 
the mass barrier, the dashed line the effect of a CP violating condensate field on 
the bubble wall. For the case shown, the net barrier is monotonic and the WKB 
reflection of particles shows no CP violation. 

Figure 2: As in Figure 1, but for the case where the net barrier is non-monotonic and the 
WKB reflection of particles does show CP violation. 

Figure 3: 'Decay' processes relevant for the propagation of a chiral tau lepton asymmetry 
through the plasma. 

Figure 4: 'Decay' processes relevant for the propagation of a chiral top quark asymmetry 
through the plasma. 

Figure 5: Vector boson t-channel exchange diagrams relevant for computation of the diffu- 
sion constants of different fermion species. 
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